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ABSTRACT 


The  Local  Emergency  Operating  System  (LEMOS)  is  a multi-program  model  within 
a system  of  models  which  are  designed  as  a part  of  a computer  based  system  to 
evaluate  local  operating  systems.  The  Defense  Civil  Preparedness  Agency 
Computation  Center  is  a co-developer  with  RTI  in  the  development  of  ADS/LEMOS. 
Effort  by  RTI  during  the  past  year  has  centered  around  the  development  of  the 
control  and  transportation  submodels  and  at  the  same  time  has  continued  to  improve 
the  procedures  within  other  submodels. 

This  report  describes  the  essential  features  of  the  control  and 
transportation  models.  In  addition,  a brief  discussion  is  included  covering  the 
application  of  the  ADS/LEMOS  model  to  the  local  CD  planning. 

The  report  concludes  that  emphasis  should  shift  now  from  design  to 
development  but  with  the  specific  objective  of  gaining  support  for  continued 
evolution  to  an  operational  state. 

Adequate  developmental  testing  is  recommended  before  any  demonstrations  of 
the  planning  role  are  undertaken. 
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DETACHABLE  SUMMARY 


The  Local  Emergency  Operating  System  (LEMOS)  is  a series  of  interrelated 
computer  programs  which  operate  as  a part  of  a larger  system  of  programs  and  manual 
procedures  designed  to  "Test  and  Evaluate  Local  Operating  Systems"  (TELOS). 

Simulation  of  local  operations  is  practicable  by  the  use  of  computer-based 
models  provided  the  user  avoids  too  much  detail  which  renders  the  simulation 
time-consuming,  and,  therefore,  expensive.  The  authors  believe  that,  while  the 
current  TELOS  design  is  seemingly  complex,  it  has,  in  fact,  avoided  this  pitfall 
and  achieved  a suitable  balance  between  too  little  and  too  much  detail.  The  level 
of  detail  adopted  is  thought  to  give  enough  realism  or  believability  and,  yet,  not 
require  endless  data  processing.  This  report  is  the  latest  in  a series  of  reports 
describing  the  evolution  of  LEMOS.  The  most  recent  effort,  culminating  in  this 
report,  includes  the  control  interface  development  and  the  transportation  model 
development  described  in  Section  II  and  III,  respectively.  While  the  interface 
between  ADS  and  LEMOS  is  believed  to  be  satisfactory,  it  has  not  been  tested  to 
verify  its  performance. 

Command  and  control  over  local  civil  defense  operations  would  be  achieved,  in 
the  real  world,  through  an  Information  system  utilizing  communication  networks 
connecting  command  and  control  centers.  In  the  LEMOS  model  command  and  control  is 
achieved  through  the  definition  of  policies,  priorities,  and  prohibitions  in  coded 
form  within  the  control  file.  The  absence  of  a communication  submodel  as  is  the 
present  case  is  tantamount  to  the  assumption  that  communications  are  "perfect". 
Since  the  prototype  LEMOS  model  does  not  Include  communications,  the  command  and 
control  model  is  defined  as  a control  file  (containing  discrete  values  for  selected 
policies,  priorities,  and  prohibitions)  and  the  Generalized  Executive  Control 
(GENEC)  which  activates  the  scenario. 

The  GENEC  system  enables  a tape  of  executable  modules  (load  and  go  modules  of 
the  major  sub-programs  in  the  LEMOS  system  plus  other  modules,  as  required)  to  be 
executed  repeatedly  in  any  sequence  desired  and  controlled  by  a second  parameter 
tape.  The  system  specifically  Iterates  between  the  countermeasures  and  damage 
assessment  models,  continuing  through  a number  of  time  periods  under  control  of 
GENEC.  Thus,  the  basic  scenario  may  be  embodied  in  GENEC  by  the  controls  residing 
In  the  Control  File  accessed  through  it.  Manual  evaluation  of  system  outputs  is 
contemplated  at  this  time.  A mechanized  output  data  processor  may  be  needed  to 
analyze  the  relatively  large  output  from  both  models  with  respect  to  the  particular 


role  of  TELOS.  On  the  basis  of  this  evaluation,  controls  for  further  scenarios  may 
be  determined  and  implemented  through  the  use  of  GENEC.  Throughout  the  course  of 
the  simulation,  reports  may  be  generated  to  measure  the  status  of  the  system  under 
test  as  an  aid  in  the  evaluation  of  local  emergency  operations. 

The  major  effort  during  this  contract  period  has  been  to  develop  the 
transportation  submodel  or,  more  accurately,  the  "shortest  path  algorithm".  This 
submodel  is  described  in  considerable  detail  and  is  viewed  as  a significant 
addition  to  LEMOS. 

If  LEMOS  is  to  have  a practical  role  in  local  CD  planning,  a pi  an  by  which 
this  role  can  be  attained  is  essential  at  this  time.  It  may  be  fairly  stated  that 
the  Research  phase  of  RDT  and  E is  nearing  completion.  The  Development,  Test,  and 
Evaluation  phases  should  begin  immediately.  During  these  phases,  selected  Federal, 
State,  and  Local  planners  should  be  enlisted  to  participate  in  evolving  and 
conducting  them.  It  is  absolutely  imperative  that  they  have  an  important  impact  on 
the  final  configuration  of  the  simulation  system  before  it  is  used  operationally  to 
improve  the  local  CD  planning  function. 

DCPA  is  strongly  urged  to  develop  appropriate  multi-year  program  plans  to 
conduct  Development,  Test,  and  Evaluation  phases  of  the  TELOS  system  with 
participation  by  Federal,  State,  and  Local  planning  personnel.  Special  emphasis 
should  be  placed  during  DT  and  E phases  on  the  use  of  Case  Study  Areas.  These 
plans  should  include  the  addition  of  the  following  elements  during  the  Development 
phase: 

Local  plans  Pre-processor, 

Interactive  Control  Procedures, 

Fire  Spread  Model  (developed  by  not  interfaced  with  ADS), 

Special  Resource  Damage  Submodel, 

Communications  Submodel, 

Utilities  Network  Submodel, 

Medical/Epidemiology  Submodel. 


LOCAL  EMERGENCY  OPERATIONS  SYTEM  (LEMOS) 


I.  INTRODUCTION 


A.  General 

The  Local  Emergency  Operating  System  (LEMOS)  is  a series  of  interrelated 
computer  programs  which  operate  as  a part  of  a larger  system  of  programs  and 
manual  procedures  designed  to  Test  and  Evaluate  Local  Operating  Systems  (TELOS). 
Figure  1 portrays  the  general  configuration  of  the  larger  system.  It  can  be  seen 
that  the  executive  control,  damage  assessment,  and  countermeasure  segments  of  that 
system  have  a particularly  close  relationship  to  each  other.  This  report  is  the 
latest  in  a series  of  reports  describing  the  evolution  of  LEMOS.  The  most  recent 
effort,  culminating  in  this  report,  includes  the  control  interface  development 
described  in  Section  II  and  the  transportation  model  development  described  in 
Section  II  and  III,  respectively.  While  the  interface  between  ADS  and  LEMOS  is 
believed  to  be  satisfactory,  it  has  not  been  tested  to  verify  its  performance. 
Future  work  should  concentrate  on  obtaining  successful  performance  tests  between 
the  various  segments  identified  in  Figure  1 and  on  performing  case  studies 
employing  TELOS  as  a planning  tool. 

First,  in  order  to  use  TELOS  as  a planning  tool,  appropriate  procedures  must 
be  developed  to  translate  existing  planning  documents  into  machine  readable  inputs 
and  to  convert  computer  outputs  into  suitable  planning  documents.  An  approach  to 
this  problem  is  described  briefly  in  the  subsection  entitled  "Local  CD  Planning." 

Second,  assuming  that  local  planning  methodology  can  be  made  compatible  with 
the  TELOS  concept,  the  local  planners  must  understand  the  simulation 
characteristics  of  the  model  and  know  how  to  use  it  to  produce  probable  outcomes 
that  they  believe  represent  their  real-world  situation.  This  methodology  Is 
discussed  under  subsection  C. 

And  third,  the  local  planners  need  to  learn  to  interpret  and  evaluate 
simulation  outcomes  in  terms  that  allow  them  to  make  decisions  regarding  the 


Figure  1.  Overall  TELOS  Description 
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status  of  local  Civil  Defense  represented  by  the  simulation  to  achieve  the  most 
effective  use  of  available  or  anticipated  resources.  The  last  subsection  in  this 
section  is  addressed  to  the  evaluation  function.  However,  it  does  not  presume  to 
cover  this  area  adequately  as,  the  authors  believe,  this  function  must  be 
developed  in  close  cooperation  with  local  officials.  For  without  their  input  and 
sanction,  TELOS  will  not  achieve  success  as  a planning  tool. 

B . Local  CO  Planning 

Planning  is  the  primary  responsibility  of  management.  In  fact,  the  local  CD 
Director  must  be  considered  to  be  the  chief  planner  and  must  devote  a very 
significant  part  of  his  time  to  discharging  this  responsibility.  Not  only  should 
he  have  written  plans  which  give  qualitative  evidence  that  he  is  fulfilling  this 
duty,  but,  they  must  be  susceptable  to  quantitative  evaluation  to  provide  adequate 
assurances  that  they  represent  the  best  available  plans.  In  short,  higher 
authorities  need  objective  methods  for  the  evaluation  of  local  CD  planning. 

In  the  era  of  detente'  and  SALT  agreements,  the  most  cost-effective  CD 
activity  that  can  be  employed  is  planning.  It  is  the  least  costly  measure  with 
the  highest  possible  payoff  and  should  precede  all  other  preparedness  activities. 

Current  Federal  guidance  provides  direct  assistance  to  the  local  planner; 
however,  it  demands  little  from  him  to  give  the  necessary  assurances  that  local 
planning  will  make  the  best  use  of  available  resources.  Nor  does  it  suggest 
reasonable  utilization  of  additional  resources  should  they  become  available. 

The  TELOS  system  of  computer  programs  offers  a means  for  measuring,  through 
simulation,  the  ability  of  local  plans  to  meet  nuclear  disaster  scenarios  deemed 
probable  for  those  areas  to  which  they  apply. 

At  present,  TELOS  has  not  been  finalized  or  tested  in  the  planning  role. 
Furthermore,  it  should  not  be  completed  until  local  planners  can  participate  in 
finalizing  the  interface  procedures  and  pass  judgment  on  its  merits  as  a planning 


tool.  TELOS  must  not  be  employed  to  denegate  the  adequacy  of  the  current  level  of 
local  planning.  In  general,  most  planners  are  painfully  aware  of  the  shortcomings 
of  their  plans.  It  should  be  used  to  stimulate  better  planning  by  revealing 
operational  inadequacies  and  reporting  on  resource  needs.  The  basic  method 
adopted  in  TELOS  is  one  of  simulating  local  operations  with  finite  resources. 

C.  Simulation  of  Local  Operations 

Simulation  of  local  operations  is  practicable  by  the  use  of  computer-based 
models  provided  the  user  avoids  too  much  detail  which  renders  the  simulation 
time-consuming,  and,  therefore,  expensive.  The  authors  believe  that,  while  the 
current  TELOS  design  is  seemingly  complex,  it  has,  in  fact,  avoided  this  pitfall 
and  achieved  a suitable  balance  between  too  little  and  too  much  detail.  The 
transportation  submodel  described  in  Section  III  typifies  this  belief.  The 
network  could  have  been  more  or  less  detailed.  The  level  of  detail  adopted  is 
thought  to  give  enough  realism  or  believability  and,  yet,  not  require  endless  data 
processing. 

The  development  of  a believable  simulation  will  depend  on  the  effective 
integration  of  the  local  planner's  ideas  into  the  simulation  process.  They  are 
the  ones  who  must  adopt  the  concept,  if  TELOS  is  to  succeed  as  a planning  tool, 
and  they  are  less  likely  to  be  advocates  of  the  system,  if  they  are  not 
participants  in  its  final  development.  The  TELOS  submodels  developed  to  date 
represent  the  main  frame  of  the  potential  structure.  The  elements  yet  to  be 
developed  are  the  elements  upon  which  the  planner  will  have  the  highest  impact. 
These  elements  may  be  loosely  characterized  as  the  input  and  output  data 
processors.  Assuming  that  appropriate  input  and  output  interfaces  are  developed 
with  local  planners  and  their  acceptance  of  the  system  attained,  planning  through 
simulation  can  proceed. 


Initially,  a local  area  may  be  described  by  its  physical  resources,  its 
countermeasure  plans,  and  the  various  probabilistic  attack  scenarios.  Simulation 
runs  may  be  made  to  determine  outcomes.  Information  gained  during  these  runs  may 
suggest  improvements  in  the  countermeasure  plans.  Revised  plans  may  be  rerun  to 
verify  the  predictions.  Quantitive  measures  gained  during  the  process  provide 
means  for  optimumizing  local  losses  and  resource  allocations. 

The  models  are  sufficiently  flexible  to  accommodate  almost  any  local 
configuration  and  may  be  run  as  may  times  as  time  and  effort  permit  to  achieve  a 
balanced  set  of  plans  uniquely  applicable  to  that  area. 

D.  Planning  Evaluation 

Since  each  local  CD  planner  is  responsible  for  his  plans,  TELOS  cannot  assume 
the  evaluative  role  for  him.  Therefore,  the  outcome  interface  must  contain 
several  alternative  means  for  evaluating  simulation  outcomes  as  they  may  be 
influenced  by  local  plans.  The  planner  should  have  not  only  the  choice  from  among 
these  alternatives  but  should  be  able  to  "tune"  the  selected  method  to  meet  his 
objectives  wherever  practicable.  Measures  for  normalizing  outcome  measures  is 
especially  important. 

Previous  reports  have  described  readiness  and  benefit  measures  which  TELOS 
may  produce  as  a result  of  the  simulation  of  local  operations.  Additional 
measures  may  be  generated  if  the  local  planner  considers  them  essential  to  the 
proper  evaluation  of  his  plans. 

Two  elements  of  local  planning  are  considered  in  the  next  two  sections  of 
this  report.  They  are  important  considerations  in  any  local  planning  simulation 
and  are  described  as  a part  of  the  development  of  the  ADS/IEMOS  subsystem. 


5 


II.  COMMANO  AND  CONTROL  MODEL 


A.  Introduction 

Command  and  control  over  local  civil  defense  operations  would  be  achieved,  in 
the  real  world,  through  an  information  system  utilizing  communication  networks 
connecting  command  and  control  centers.  In  the  LEMOS  model  with  a communication 
(COM)  submodel,  command  and  control  is  achieved  through  the  definition  of 
policies,  priorities,  and  prohibitions  in  coded  form  within  the  control  file.  The 
absence  of  a communication  submodel  is  tantamount  to  the  assumption  that 
communications  are  "perfect".  In  this  view,  the  presence  of  a COM  submodel 
permits  the  evaluation  of  a delay  between  an  event  and  the  generation  of  data 
(information),  between  its  generation  (transmission)  and  reception,  or  between  its 
reception  and  decision-making.  Loss  of  information  is  equivalent  to  an  infinite 
del  ay. 

Since  the  prototype  LEMOS  model  does  not  include  communications,  the  command 
and  contol  model  is  defined  as  a control  file  (containing  discrete  values  for 
selected  policies,  priorities,  and  prohibitions)  and  the  Generalized  Executive 
Control  (GENEC)  which  activates  the  scenario. 

8 . Generalized  Executive  Control 

The  Generalized  Executive  Control  (GENEC)  system  represented  by  Figure  2 was 
developed  by  the  Defense  Civil  Preparedness  Agency  Computer  Center  (DCPACC)  and 
was  designed  to  run  on  the  COC  3600  at  DCPACC  under  the  SCOPE^  operating  system. 
The  CDC  3600  and  SCOPE  system  is  a single  program  execution  computer  system  with 
no  inherent  provisions  for  sequentially  repeating  a set  of  programs.  The  ADS  and 
LEMOS  series  of  programs  should  be  executed  as  a system  and  not  as  separate  runs 
due  to  the  interactive  and  iterative  nature  of  the  process.  The  GENEC  system 
was  specifically  developed  to  accomplish  this  task.  Its  operation  is  shown 
di agrammatical ly  in  Figure  2 and  in  an  alternate  system  concept,  would  be 

* software  system  operating  on  the  CDC-3600. 
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Program*  In  R*  locacsb  La  Form 

(LOAD- AND- GO)  Paramacer*  for  Running  Program* 


Contain: 

Numbar  of  Timas 
Numbar  of  Jobs  (Program) 

Position  of  First  Job  to  b*  Run  on  LGO  Tapa 
Position  of  First  Paramatar  Sat  on  PAR  Tap* 


Figure  2.  Conceptual  Model  of  GENEC 

interactive  and  capable  of  being  dynamically  modified  as  it  executes  successive 
passes. 

GENEC  is  initiated  by  control  card  input,  which  tells  GENEC  the  number  of 
times  that  the  job  stream  is  to  be  repeated,  the  number  of  jobs  in  the  stream,  the 
position  of  the  first  job  on  the  program  tape  and  the  position  of  the  first  set  of 
parameters  on  the  parameters  tape.  In  addition,  GENEC  will  be  able  to  pass 
approximately  five  words  of  data  from  one  routine  to  the  next  routine  executed  by 
the  memory  mapping  technique. 
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C.  ADS/LEMOS/GENEC  Operation 

The  overall  flow  of  the  TELOS  system  is  under  the  control  of  manual 
intervention  and  GENEC  as  outlined  in  Figure  2. 

The  GENEC  system  enables  a tape  of  executable  modules  (load  and  go  modules  of 
the  major  sub-programs  in  the  LEMOS  system  plus  other  modules,  as  required,  from 
AOS)  to  be  executed  in  any  sequence  and  executed  repeatedly,  if  desired,  as 
controlled  by  a second  control  or  parameter  tape. 

The  system  specifically  iterates  between  the  countermeasures  and  damage 
assessment  models,  continuing  through  a number  of  time  periods  under  control  of 
GENEC.  Thus,  the  basic  scenario  may  be  embodied  by  GENEC  and  controls  residing 
in  it  and  in  the  Control  File  accessed  through  it.  Manual  evaluation  of  system 
outputs  is  contemplated  at  this  time.  A mechanized  output  data  processor  may  be 
needed  to  analyze  the  relatively  large  output  from  both  models  with  respect  to  the 
input  and  the  particular  role  of  TELOS.  On  the  basis  of  this  evaluation,  controls 
for  further  iterations  may  be  determined  and  implemented  through  the  use  of  GENEC, 
and  a new  cycle  begins.  Throughout  the  course  of  the  simulation,  reports  may  be 
generated  to  measure  the  status  of  the  system  under  test  as  an  aid  in  the 
evaluation  of  local  emergency  operations. 

The  main  data  linkage  between  all  of  the  various  LEMOS  submodels,  as  shown  in 
Figure  3,  is  the  Control  File,  although  some  13  file  types  are  used  together  or 
individually  as  transitions  between  programs.  The  Control  File  is,  however,  only 
modified  by  the  manual  intervention  exterbak  to  ADS/LEMOS. 

0.  GENEC/Control  File  Interaction 

As  a part  of  each  sub-program  in  the  LEMOS  system,  there  is  a set  of 
parameters  that  determine  the  type  of  run,  type  of  input,  and  type  of  output  that 
a module  is  to  provide. 

The  "essential"  parameters  are  used  to  control  the  functions  of  each  of  the 
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(Being  Drafted) 


Figure  3.  Local  Operating  System  Overview 


I 
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programs  in  LEMOS.  There  are  three  "essential"  parameters  which  are  input  to  each 
program,  called  RUN-SWITCH,  TST-SWITCH,  and  PRNT-SWITCH.  The  first  of  these, 
RUN-SWITCH,  is  generally  used  to  select  which  major  functions  a program  should 
perform.  The  second  parameter,  TST-SWITCH,  is  generally  used  for  aid  in  debugging 
or  testing  a program.  The  last  parameter,  PRNT-SWITCH,  is  generally  used  to 
select  printout  options  within  a program,  e.g.,  which  table(s)  to  list. 

A GENEC  modification  provides  a way  to  pass  a few  parameters  from  one  module 
to  another  (five  48  bit  words)  and,  therefore,  the  execution  of  the  full  LEMOS 
system  can  be  previous  history  dependent  as  well  as  having  its  basic  execution 
being  dependent  on  the  planned  scenario. 

The  control  file  and  procedures  to  access  it  are  contained  in  a copy  library 
which  is  invoked  at  the  time  of  program  compilation.  This  procedure  assures 
uniformity  in  the  application  of  policies,  priorities,  and  prohibitions  to  all 
procedures. 

All  programs  in  the  LEMOS  model,  including  the  transportation  model  described 
in  the  next  section,  have  access  to  these  control  procedures. 


III.  TRANSPORTATION  MODEL 


A.  General 

The  major  effort  during  this  contract  period  has  been  to  develop  the 
transportation  submodel  or,  more  accurately,  the  "shortest  path  algorithm".  This 
submodel  is  described  in  considerable  detail  in  Appendix  B and  discussed  in  the 
following  sub-sections.  It  is  viewed  as  a significant  addition  to  LEMOS  and, 
indeed,  to  the  entire  ADS/LEMOS  system. 

8.  Network  Description 

Freeways  and  major  arteries  are  the  basis  of  the  current  development  of  the 
transportation  network.  Residential  and  feeder  streets  are  not  included  as  they 
unnecessarily  expand  the  network  to  unmanageable  proportions.  Intersections  of 
highways  and  arteries  define  nodes  in  the  network.  Additional  nodes  could  be 
added,  if  desired,  at  non-intersection  points.  The  segment  between  two  nodes  is  a 
link  in  the  network.  Nodes  are  identified  by  unique  numbers.  Connectivity  of 
Unit  Areas  is  determined  by  the  occurrence  of  equal  node  numbers.  Links  are 
identified  by  numbers  and  names.  Link  identification  is  not  involved 
computationally  in  the  current  computerized  transportation  sub-model;  it  is  used 
only  to  aid  in  creating  the  links  file  data  base.  The  transportation  network  is 
described  by  the  Links  File.  The  data  in  this  file  is  used  to  create  the  matrix 
of  distances  between  directly  linked  nodes  in  a unit  area  network.  Each  record  of 
this  file  defines  a link,  primarily,  in  terms  of  its  two  terminal  nodes,  the 
distance  in  miles  between  the  nodes,  and  the  type  of  link  (i.e.,  one-way  or  both 
ways).  The  two  nodes  are  identified  as  a backward  node  and  a forward  node.  The 
distinction  is  not  critical  except  in  the  case  of  one-way  streets  and  when  the 
distance  from  one  node  to  the  other  differs  from  the  converse  in  which  two 
different  links  records  are  prepared.  The  link  type  is  used  to  distinguish 
between  one-way  and  two-way  links.  It  could  also  be  used  to  describe  other 


transportation  networks,  i.e.,  railroads  or  waterways.  However,  our  present 
effort  is  limited  to  highway  networks.  In  earlier  work,  the  forward  node  was 
defined  as  follows:  if  the  angle  measured  from  due  North  taken  about  one  of  the 

two  nodes  to  the  straight  line  connecting  the  nodes  is  greater  than  180  degrees, 
that  node  is  the  forward  node.  From  this  it  follows  that  the  other  node  is  the 
backward  node  and,  if  the  angle  is  less  than  180  degrees,  the  converse  holds. 


illustrates  the  definition  of  forward  and  backward  nodes.  The  indicated  angle,  a, 
about  node  A is  less  than  180  degrees,  therefore,  A is  the  backward  node  and  B is 
the  forward  node. 

In  the  current  version,  the  definition  is  retained,  but  the  emphasis  is  on 
the  distance  and  direction  of  travel  between  the  two  nodes  and  the  associated  link 
codes.  The  table  in  Appendix  B,  Input-Output  Description  Section,  summarizes  the 
codes  used  in  the  links  file. 

In  the  present  version  of  the  computerized  Transportation  Submodel,  no  test 
is  made  on  the  first  character  of  the  link  code.  Thus,  only  major  arteries  (i.e., 
roadways)  should  be  included  in  the  Links  File.  If  other  types  are  included,  it 

should  be  borne  in  mind  that  a constant  speed  of  45  mph  is  used  to  convert 

distances  to  travel  time,  as  coded  in  the  Links  File  and  adjusted  by  the  Basic 

Operating  Situation  (BOS)  in  the  Resource  File. 

The  table  in  Appendix  B,  Input-Output  Description,  describes  the  content  of 
the  records  in  the  Links  File. 


In  the  previous  work  Ref  1 , four  levels  of  networks  coinciding  with  the 
organizational  structure  were  devised  because  the  estimated  number  of  links  in  the 
anticipated  networks  exceeded  computer  capacity  when  the  standard  procedure  of 
storing  the  distance  matrix  in  computer  core  was  used.  The  current  computerized 
version,  using  out-of-core  storage  of  the  matrix,  can  deal  with  relatively  large 
networks.  The  reasons  for  retaining  the  network-level  and  node-level  codes  in  the 
Links  File  data  are  to  provide  visual  checks  in  data  preparation  and  to  provide 
the  means  necessary  to  view  the  transportation  problem  as  sector  to  sector  instead 
of  unit  area  to  unit  area,  if  this  should  ever  become  desirable  in  the  future. 
Essentially,  the  only  changes  necessary  would  be  to  remove  the  '999'  divider  cards 
in  the  Links  File  between  those  unit  areas  to  be  agregated  and  to  change  the 
coding  in  tests  of  the  network  identification  number  to  test  the  value  of  the 
variable  corresponding  to  the  appropriate  level.  Linkage  between  higher  level 
networks  would  be  determined  by  equality  of  node  identification  numbers  having  a 
node-level  code  appropriate  for  the  network  level. 

Node-level  codes  and  their  meanings  are  summrized  as  follows: 

Node  Level 

Code  Meaning 

1 Interior  to  unit  area;  not  shared  with  other  areas 

2 Shared  between  unit  areas 

3 Shared  between  sectors 

4 Shared  between  groups 

5 Shared  between  EOC's 

6 On  the  boundary  of  zones. 

A node  with  a level  code  of  4 is  shared  by  (at  least)  two  groups,  i.e.,  it  is 
on  the  boundary  between  two  groups.  The  node  may  occur  as  a boundary  node 
between  some  sectors,  constituents  of  those  groups,  and  some  unit  areas. 
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constituents  of  those  sectors.  The  value  assigned  the  node-level  code  should 
correspond  to  the  highest  level  occupied  by  the  node.  This,  however,  is  not  a 
requirement  for  the  correct  execution  of  the  current  version  of  the  computer 
program. 

C.  Minimum  Travel  Time  (Path)  Between  Unit  Areas 

The  computerized  Transportation  Submodel  computes  the  paths  of  minimum  travel 
time  between  all  unit  areas  and  creates  the  TVL-REC  file.  The  problem  is  defined 
and  solved  in  terms  of  distances,  and  the  conversion  from  distance  to  travel  time 
is  made  by  assuming  a speed  of  45  miles  per  hour.  The  input  data  to  the  computer 
model  consist  of  two  files,  the  Links  File  describing  each  unit  area  network  in 
its  initial  state,  and  a Control  File,  which  identifies  the  numbers  of  the  unit 
area  networks  selected  for  processing.  The  computer  program  has  been  written  to 
accommodate  unit  area  networks  containing  a maximum  of  75  nodes  each  and  a network 
comprising  a maximum  of  400  unit  areas.  It  has  been  tested  on  a network  of  8 
unit  areas  containing  from  6 to  26  nodes. 

There  are  two  major  steps  in  the  model.  The  first  is  a "one-for-all"  step 
that  computes  shortest  paths  between  all  nodes  in  a unit  area  network  and  the 
distances  between  directly  linked  unit  area  networks.  This  step  used  the  Link  and 
Control  Files  previously  mentioned.  The  second  step  uses  the  Basic  Operating 
Situation  of  each  unit  area  to  alter  these  distances  and,  then,  to  compute  the 
shortest  paths  between  all  unit  areas.  The  second  step  can  be  repeated  as  often 
as  necessary  as  BOS  changes  occur  over  time  in  a given  scenario. 

Both  steps  use  Floyd's  Algorithm  to  compute  the  shortest  paths  between  all 
nodes.  The  major  difference  between  the  two  steps  in  the  application  of  the 
algorithm  Is  that  in  the  first  step  (distance  between  nodes  in  a unit  area)  the 
entire  distance  matrix  is  stored  in  computer  core  memory  while  the  second  step 
takes  advantage  of  the  fact  that  the  algorithm  operates  upon  only  two  rows  of  the 
distance  matrix  at  a time.  Thus,  only  two  rows  of  the  distance  matrix  need  be  in 
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core  at  any  time.  The  remainder  of  the  matrix  is  stored  on  tape  or  disk. 
Substitution  of  tape  or  disk  for  computer  core  memory  permits  the  direct 
processing  of  extremely  large  matrices,  i.e.,  large  numbers  of  unit  areas,  which, 
if  stored  in  core,  would  exceed  the  capacity  of  the  computer  to  be  used.  This 
ability  to  compute  directly  the  shortest  paths  between  unit  areas  obviates 
resorting  to  the  previously  described  hierarchical  scheme  in  which  unit  areas  were 
aggregated  into  sectors,  sectors  into  groups,  and  so  on,  where  shortest  paths  were 
computed  between  the  constituents  of  each  hierarchical  level. 

Two  unit  area  networks  are  directly-linked  if  they  share  one  or  more  nodes  in 
common.  Such  nodes  are  called  boundary  nodes.  The  distance  between  two  directly- 
linked  unit  area  networks  from  A to  B is  defined  as  the  average  distance  taken 
over  shortest  paths  from  all  nodes  in  A to  all  nodes  in  B via  the  boundary  nodes 
between  A and  B.  If  A and  B are  not  directly-linked,  the  distance  between  them  is 
set  at  a large  number  for  computational  purposes.  The  distance  from  A to  B is  not 
necessarily  equal  to  the  distance  from  B to  A because  of,  for  example,  the 
possible  occurrence  of  one-way  streets.  In  application,  two  distance  matrices  are 
maintained.  In  the  first  matrix,  each  element  represents  the  average  distance 
(over  shortest  paths)  from  all  nodes  in  a given  unit  area  to  the  boundary  nodes  it 
shares  with  some  directly-linked  unit  area.  The  elements  of  the  second  matrix  are 
the  average  distances  from  the  common  boundary  nodes  to  all  other  nodes  in  the 
other  unit  area.  Each  element  of  the  final  matrix  of  distances  between  directly- 
linked  unit  areas  is  the  sum  of  the  corresponding  elements  of  these  two  matrices 
divided  by  the  corresponding  number  of  common  boundary  nodes. 

In  a time  scenario,  damage  is  inflicted  upon  (or  removed  from)  a unit  area 
resulting  in  a change  in  travel  time  within  the  area.  The  status  of  the  unit  area 
is  not  determined  to  any  finer  resolution  down  to  the  level  of  individual  blocks. 
Thus,  the  damage  is  assumed  to  be  distributed  uniformly  over  the  unit  area.  All 
initial  state  shortest  paths  within  a unit  area,  in  effect,  are  increased  by  a 
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factor  greater  than  1,  if  the  state  of  the  area  has  deteriorated,  or  decreased  by 
a factor  less  than  1,  if  conditions  have  improved.  This  apparent  change  in 
distance  results  in  an  observed  change  in  travel  time. 


Values  of  BOS  range  from  1 to  9 and  indicate  the  status  of  the  environment  in 
a unit  area.  The  value  includes  the  effects  of  debris,  fires,  and  radiation. 
Somewhat  arbitrarily,  the  increase  in  travel  times  associated  with  an  increase  in 
BOS  is  functionally  represented  by 


f = 2 


(BOS-1) 


Values  of  BOS,  their  meaning,  and  the  associated  increase  in  travel  times  are 


summarized  as  follows: 


Meaning* 


tadio- 


Values  Debris  Fires  activity  Increase  in  Travel  Time 


* N = Negligible  M = Moderate  S = Severe 


X 128 
X 256 


The  intent  is  to  increase  the  travel  time  in  an  area  with  a BOS  of  9 so  that 
the  shortest  paths  algorithm  is  unlikely  to  select  that  area  as  an  intermediate 
node.  In  other  words,  such  areas  are  to  be  avoided.  Alternative  weighting 
systems  may  be  used  if  the  results  warrant  changing  the  method  adopted  herein. 
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While  it  is  not  necessary  to  recompute  the  shortest  paths  within  a unit  area, 


I 
I 

f since  all  distances  are  multiplied  by  the  same  factor,  it  is  necessary  to 

recompute  the  shortest  paths  between  unit  areas.  To  illustrate,  consider  the  case 
where  a single  unit  area,  Q,  has  sustained  damage,  i.e.  an  increase  in  BOS,  in  a 
network  comprising  a total  of  n unit  areas.  Basically,  two  situations  can  occur. 

The  first  situation  exists  where  Q is  the  origin  or  the  destination  in  the 
shortest  path  between  two  unit  areas.  These  shortest  paths,  containing 
intermediate  unit  areas  which  have  not  sustained  a change  in  BOS,  are  only  changed 
by  an  additive  amount  resulting  from  the  increase  sustained  by  Q,  which  is  at  the 
end  of  the  paths.  Thus,  it  is  not  necessary  to  employ  an  algorithm  to  resolve  the 
shortest  paths  problem  in  this  situation. 

The  second  situation  involves  Q as  an  intermediate  area  in  the  shortest  path 
between  two  other  areas.  In  this  situation,  a shorter  path  may  result  by  using 
some  area  other  than  Q as  an  intermediate  area.  Thus,  it  is  necessary  to  resolve 
the  shortest  paths  problem  in  this  situation.  A description  of  the  various 
techniques  considered  for  application  in  this  situation  follows. 

In  addition  to  Floyd's  algorithm  (which  solves  for  all  shortest  paths  from 
all  nodes  in  a network  to  either  one  or  all  other  nodes),  another  efficient  method 
is  Dijkstra's  algorithm,  which  solves  for  either  the  shortest  path  between  two 
specified  nodes  or  the  shortest  paths  from  a specified  origin  to  all  destinations. 
Identification  of  those  shortest  paths,  where  Q is  an  intermediate  node,  would 
then  provide  a list  of  origin/destination-specific  paths  to  be  resolved  using 
Dijkstra's  algorithm.  The  disadvantage  to  this  approach  lies  in  the  storage  and 
retrieval  by  rows  of  not  only  the  matrix  of  shortest  paths  but  also  the  policy 
matrix  used  initially  by  Floyd’s  algorithm  to  solve  for  all  shortest  paths  between 
all  nodes.  Essentially,  the  disadvantage  is  that  an  excessive  amount  of  computer 
input/output  time  would  be  consumed  in  reading  from  storage  various  rows  of  these 
two  matrices  in  a random,  non-predictable  order.  The  basis  for  this  conclusion 
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will  be  discussed  in  greater  detail.  However,  as  a consequence  of  it,  the  problem 
can  be  resolved  by  a re-application  of  Floyd's  algorithm  in  approximately  the  same 
amount  of  computing  time  required  for  an  application  of  Dijkstra's  procedure  with 
the  added  benefit  of  reduced  program  coding,  size,  and  maintenance. 

Designate  the  matrix  of  distances  between  n nodes  by 


D = (d  ) , (i,j  = 1,2, ,n) 

i.j 

In  particular  the  matrix  of  distances  between  directly-linked  nodes 
designated 

0 0 


0 

where:  d 

0 


D = (d  ) 
i.j 

°°  if  nodes  i and  j are  not  directly-linked 


and  d < °°  if  nodes  i and  j are  directly-linked, 

i J 


is 


The  policy  matrix  is  denoted  by 

P = (P  ) . (i.j  = 1.2, ,n) 


i.j 

and,  in  particular,  the  initial  policy  matrix  is 


0 0 
P = (p  ) 


where:  p = 0 if  i and  j are  not  directly-1 inked 

i.j 

0 

and  p = j if  i and  j are  directly- linked. 

i.j 


In  Floyd's  algorithm  each  node  k = l,2,...,n  is  utilized  in  turn  as  an 
intermediate  node  in  the  path  between  all  other  nodes  i.j  * l,2,...n.  This  means 
that  the  trial  distance  using  node  k is  computed  using 


i+k  k k 

d = d + d 
i,j  1,k  k,j 


If  the  trial  distance  is  less  than  the  previous  distance  d , then  the  trial 

i .3 
k 

distance  replaces  d in  the  distance  matrix  D,  i.e., 
i,J 

k k k 

d + d d 

i.k  k , j i , j 


and,  also 


P + 9 
i,k  i,j 


after  n interations  (k  = 0,1,2,  .f . ,n-l)  on  the  matrix  D,  the  final  distance  matrix 


of  shortest  paths. 


n n 

0 = (d  ) , 

1 .J 


and  the  final  policy  matrix. 


n n 

P = (P  ) , 

i,j 


are  obtained. 


Note  that  in  the  computation  of  the  trial  distances, 

k+1  k k 

d = d + d , 

1.J  l.j  k,  j 


and  the  comparisons  to  d , only  elements  of  two  rows,  the  ith  and  the  kth, 

1,j 

are  involved.  It  is  this  feature  which  has  been  programmed  to  permit  the  direct 
solution  of  the  shortest  paths  between  all  unit  areas. 

The  interpretation  of  a solution  policy  matrix,  P,  (dropping  the  superscript 
notation)  is  as  follows.  For  illustration,  it  is  true  that  the  shortest  path 
between  two  nodes,  i and  j,  involves  the  following  intermediate  nodes  in  order. 


• 


i.e.,  the  policy  to  go  from  i to  m is  i to  1 to  m to  q to  r to  j,  and  also  assumes 
that  q is  a unit  area  that  will  subsequently  be  damaged.  An  abbreviated  listing 
of  the  policy  matrix  resulting  from  application  of  Floyd's  algorithm  would  be  as 
follows: 


i 

a b . . . i j . . . z 


The  partial  listing  above  of  a Policy  Matrix  is  used  to  illustrate  a 
hypothetical  policy  for  the  shortest  path  from  node  i to  node  j;  i to  1,  1 to  m, 
m to  q,  q to  r,  and  r to  j. 
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-i. 


■ v 
N- 


Examination  of  P to  obtain  the  policy,  p , would  yield  in  succession 

P = 1 . 

1 J 

P = m , 

l.j 

P ■ q , 

fn.j 

P = r , 

q. j 

and,  finally,  p 3 j . 

r, j 

Therefore,  the  shortest  path  policy  from  node  i to  node  j,  as  found  by 
examination  of  p,  is  given  by 

p = ( i , 1 ,m,q ,r , j)  . 

This  concludes  the  interpretation  of  the  policy  matrix,  P = (p  ). 

i,j 

The  following  discussions  of  the  damage  assessment  problem  will  describe  the 
computational  steps  required  to  determine  the  occurrence  of  node  q as  an 
intermediate  node;  Dijkstra's  procedure  for  determining  the  shortest  paths  between 
a specified  origin  and  a specified  destination;  an  evaluation  of  Dijkstra's 
procedure  versus  Floyd's  procedure;  and  a feature  of  Floyd's  algorithm  which  could 
be  implemented  with  possible  benefit. 


D.  The  Damage  Assessment  Problem 

Assume  that  node  q has  sustained  damage.  From  the  illustrative  example. 


it  is  true  that  in  the  solution  policy  matrix,  P,  that 


P = q 

m,r 


and  p = q . 
m,j 

That  is,  in  application,  upon  reading  the  mth  rows  of  the  policy  matrix  the  value 
of  q will  be  found  to  occur  two  times  indicating  that  q is  an  intermediate  node 
in  the  shortest  paths,  m to  r and  m to  j. 

It  is  preferable  to  solve  only  the  shortest  path,  m to  r,  since  the  policy 


P = (m,q,r) 

m,r 


occurs  in  p = (m,q,r,j)  . 


This  information  is  not  available  at  this  point;  only  p = q and  p = q are 

m,r  m,  j 

known.  To  obtain  the  full  policies,  m to  r and  m to  j,  it  is  necessary  to 

examine  the  rth  and  the  jth  columns  of  p which  is  stored  by  rows.  To  examine  p by 

columns  (or,  equivalently,  transpose  rows  and  columns),  m dual  I/O  operations 

(read  and  rewind  are  performed  on  the  complete  n x n matrix  P)  are  necessary. 

Assume  this  is  done  and  both  policies,  p(m,r)  and  ?(m,j)  are  obtained.  Recalling 

that  policies  are  stated  in  terms  of  directly-linked  nodes,  examination  of  the  qth 

0 

row  of  the  initial  policy  matrix,  p , for  those  nodes  directly-linked  to  q could 
show  that  not  only 

0 

P = r , 

q.r 

0 

but  also  p 3 j . 

qJ 
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That  is,  q is  directly-linked  to  both  r and  j and  both  shortest  paths  must  be 
solved.  If  instead 

0 

p 3 0 

q.j 

were  to  occur,  it  would  then  be  clear  to  solve  only  for  the  shortest  path  m to  r. 
Thus,  transposing  the  solution  policy  matrix  and  reading  the  initial  policy  matrix 
are  necessary  to  determine  whether  a damaged  node  is  an  intermediate  node  and  to 
determine  which  shortest  paths  should  be  recomputed. 


and 


0,  j - 1 


/ » J * 

{L  | » { 

J [ 00 , j 3 2,3, ,n  . 


The  algorithm  proceeds  by  setting  i = 1 and  comparing  the  labels  |L  } with  the  sum 

1 

of  the  kth  label  and  direct  distance  from  node  k to  node  j: 

i 

min  (L  , L + d ) -*■  L , 

j k k j 

i i,j 

for  all  j not  in  |k  l.  Examining  the  resulting  vector  of  labels  |l  } for  min  JL  } 

1 j j 

the  associated  node  index,  k , is  added  to  the  array  |k  l.  This  step  removes 

1+1  i 

that  node  from  further  consideration  and  places  it  in  the  shortest  path  from  node 
1 to  node  n.  Then,  i is  incremented 
i + 1 - 


and  the  step 

. /j  ~ 1,2, .... ,n, 

min(L,L  +d 

j 'k  k j 'j7|k( 

i 1,J  1 

is  repeated;  min  |L  t is  found;  the  associated  index  added  to  j k l;  and  so  on 
j 1 

until,  after,  at  most,  n-1  such  iterations,  min  jl  ^ is  found  to  be  the  node,  n. 

j 


F.  Evaluation  of  Dijkestra's  Versus  Floyd's  Algorithm 

In  each  of  the  i = 1,2 ,n,  m < n-1  iterations  the  direct  distance  from 

node  k to  all  other  nodes  j 3 1,2, ,n  (j  = |k  })  is  needed.  As  a result,  the 

1 i 

distance  matrix,  (d  ),  which  is  stored  as  a sequential  file  of  rows  on  disk  or 

1.J 

tape,  may  have  to  be  accessed  (rewound  and/or  partially  read  to  row  k ) as  many 

i 

as  n times  for  the  hypothetical  example  of  a single  damaged  intermediate  node. 
These  accessings  of  the  file  would  he  necessary  for  each  unique  occurrence  of  a 
node  as  a damaged  intermediate  node  and  for  all  such  nodes.  The  un-predictable 
number  of  I/O  operations  plus  the  number  required  to  transpose  the  policy  matrix 
are  the  factors  which,  in  addition  to  considerations  of  computer  programming 
simplicity  and  maintenance,  have  lead  to  the  decision  of  preferring  Dijkstra's 
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procedure  in  favor  of  Floyd's. 


G.  Possible  Benefit  from  Implementation  of  Recomputation  Procedure 

— k+1 

In  Floyd's  algorithm,  n successive  distance  matrices,  D (k=0,l,l, ,n-l) 

are  computed.  The  k+1  matrix  represents  the  shortest  paths  between  all  node  pairs 
where  each  node  i*l ,2 , . . . ,k-l,k  has  been  used  as  an  intermediate  node.  If  node  q 
is  a damaged  intermediate  node,  the  shortest  paths  may  be  a recomputed  beginning 
with  the  matrix  in  which  node  q-1  was  evaluated  as  an  intermediate  node.  In  other 
words,  all  shortest  paths  based  on  non-damaged  intermediate  nodes  1,2,..., q-1  are 
still  valid;  the  algorithm  does  not  need  to  be  repeated  for  node  indices  < q. 

The  technique  has  not  been  implemented,  but  in  principle  at  least,  the 
successive  matrices  are  saved  by  writing  them  as  a sequential  file  on  disk  or  tape 
storage.  From  a list  of  the  identification  numbers  of  unit  areas  that  have 
sustained  a change  in  BOS  since  the  last  previous  time  step  in  the  scenario,  the 
minimum  (numerical  value)  is  determined,  e.g.,  the  value  of  the  minumum  is  q.  The 
file  of  saved  distance  matrices  is  advanced  to  the  beginning  of  the  matrix 
corresponding  to  q-1.  That  distance  matrix,  then,  is  read  as  the  initial  distance 
matrix  to  begin  the  shortest  paths  algorithm  with  q as  the  first  tiral 
intermediate  node. 

There  are  several  reasons  for  not  implementing  the  technique  at  this  time. 

One  question  pertains  to  the  trade-off  between  increased  computer  I/O  time  and 
simply  recomputing  the  shortest  paths  beginning  at  node  one.  If  implemented,  the 
technique  would  require  that  the  policy  matrix  be  computed  and  saved  on  disk  or 
tape  in  addition  to  the  distance  matrix  at  each  iteration,  and  examined  if  it  is 
desired  to  begin  the  algorithm  at  the  distance  matrix  corresponding  to  the  minimum 
damaged  intermediate  node.  As  an  alternative,  the  question  of  the  node  being  an 
intermediate  node  could  be  ignored;  simply  re-apply  the  algorithm  beginning  with 
the  matrix  that  corresponds  to  the  minimum  damaged  node.  This  would  obviate  the 


necessity  of  computing/saving  the  n policy  matrices  and  transposing/examining  the 
final  pol icy  matrix. 

Another  question  originates  from  the  definition  of  distances  between 
directly-linked  unit  areas.  If  areas  A and  B are  linked  and  A is  damaged,  the 
distance  from  A to  the  boundary  between  A and  B is  changed  while  the  distance  from 
the  boundary  to  B is  not.  These  two  distances  are  stored  in  two  different  files 
and  are  added  to  obtain  the  distance  from  A to  B.  The  unresolved  questions  are  at 
what  stage  of  the  iteration  and  with  what  additional  complexity  and  increased  I/O 
time  are  the  distance  changes  and  additions  to  be  performed  and  with  what  benefit 
when  there  are  large  numbers  of  damaged  areas. 

To  derive  the  maximum  benefit  from  the  technique,  ideally,  the  area  to  be 
damaged  (say  there  is  only  one  for  illustration)  should  be  the  n-th;  only  one 
iteration  of  the  algorithm  would  be  necessary  to  resolve  the  shortest  paths.  If 
the  damaged  area  were  the  1st,  then  the  algorithm  must  be  completely  re-iterated, 
and  the  technique  has  been  of  no  benefit.  These  observations  suggest  that  the 
ordering  of  the  unit  area  identification  or  the  structuring  of  the  Links  File 
should  have  those  areas  most  likely  to  be  damaged  and  intermediate  ordered  last. 
This  criterion  has  impact  on  and  interacts  with  other  aspects  of  the  total  overall 
problem.  Since  it  could  not  be  evaluated  it  remains  an  unresolved  question. 

Finally,  data  for  8 unit  area  networks  were  available  for  checkout  and 
evaluation.  This  small  amount  of  data  would  not  permit  answering  questions  of 
feasibility  in  the  case  of  much  larger  sets  of  unit  areas. 


IV.  DISCUSSION 


The  discussion  of  the  Local  Emergency  Operating  System  (LEMOS)  is  undertaken 
in  the  context  of  its  interactions  with  AOS,  the  attack  environment,  and  fire 
spread  models.  Any  simulation  is  impossible  without  them.  Initial  discussion 
will  focus  on  the  need  for  a Development,  Test  and  Evaluation  (DT  and  E)  plan  and, 
then,  suggest  perceived  needs  for  both  ADS  and  LEMOS. 

If  LEMOS  is  to  have  a practical  role  in  local  CD  planning,  a plan  by  which 
this  role  can  be  attained  is  essential  at  this  time.  It  may  be  fairly  stated  that 
the  Research  phase  of  RDT  and  E is  nearing  completion.  The  Development,  Test,  and 
Evaluation  phases  should  begin  immediately.  During  these  phases,  selected 
Federal,  State,  and  Local  planners  should  be  enlisted  to  participate  in  evolving 
and  conducting  them.  It  is  absolutely  imperative  that  they  have  an  important 
impact  on  the  final  conf iguration  of  the  simulation  system  before  it  is  used 
operationally  to  improve  the  local  CD  planning  function. 

While  the  mainframe  of  ADS/LEMOS  is  nearly  complete,  several  members  are 
conspicuously  incomplete.  First,  ADS  is  without  a special  resources  submodel  and 

an  adequate  fire  spread  model.  Second,  LEMOS  is  without  a conmunications  module 

< 

to  approximate  realistically  the  communications  problems  within  a damaged  area. 
(The  current  model  assumes  no  communication  problems,  i.e.,  not  on  the  problem 
file  and  that  all  needed  communication  can  be  completed  as  needed.)  Third,  an 
upgraded  version  of  the  DCPA  Emergency  Medical  Model  is  needed  with  its  resource 
requirements  incorporated  into  the  LEMOS  system.  The  resulting  outputs  from  the 
epidemiology  model  are  needed  inputs  to  the  overall  system.  Last,  but  not  least, 
the  GENEC  system,  as  described  earlier,  does  not  allow  "on  line"  or  "interactive" 
modifications  to  the  run;  however,  with  the  current  advances  in  operating  systems, 
this  implementation  could  be  Incorporated  during  the  DT  and  E phases,  if  plans  are 
evolved  to  do  so. 


27 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


DCPA  is  strongly  urged  to  develop  appropriate  multi-year  program  plans  to 
conduct  Development,  Test,  and  Evaluation  phases  of  the  TELOS  system  with 
participation  by  Federal,  State,  and  Local  planning  personnel.  Special  emphasis 
should  be  placed  during  DT  and  E phases  on  the  use  of  Case  Study  Areas.  These 
plans  should  include  the  addition  of  the  following  elements  during  the  Development 
phase: 

Local  plans  Pre-processor, 

Interactive  Control  Procedures, 

Fire  Spread  Model  (developed  by  not  interfaced  with  ADS), 

Special  Resource  Damage  Submodel, 

Communications  Submodel, 

Utilities  Network  Submodel, 

Medical/Epidemiology  Submodel. 
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USER’S  GUIDE 


1.  Introduction 

This  section  provides  a general  description  of  the  procedures  required 
to  use  the  Local  Emergency  Operating  System.  It  is  assumed  that  this  usage 
will  occur  as  a segment  of  the  TELOS  model,  under  the  control  of  GENEC. 

2.  Use  of  Master  File 

The  Master  Status  File,  described  in  detail  in  Reference  1,  is  the  major 
link  between  ADS  and  LEMOS  within  TELOS.  The  Master  Status  File  contains  a 
set  of  records  for  all  resources  except  network  resources  (e.g.,  highways 
and  power  lines)  and  their  damage  states  in  all  the  unit  areas  comprising 
the  zone  being  studied.  (Unit  areas  are  defined  in  reference  1.)  These  re- 
sources include  structures,  shelter  spaces,  personnel  (both  civil  defense 
and  general  population),  and  civil  defense  assets  (including  teams,  equip- 
ment sets,  and  supplies).  Thus,  the  Master  Status  File  is  best  described 
as  a temporary  data  base  for  the  zone  being  studied. 

The  Master  Status  File  can  be  used  in  two  different  modes  by  TELOS: 
a "static"  mode  and  a "dynamic"  mode.  The  former  does  not  involve  the  use 
of  LEMOS  which  is  indicated  by  the  term  "static"  (with  respect  to  civil 
defense  countermeasures) . The  "dynamic"  mode  is  a multi-pass  run  of  the 
TELOS  using  civil  defense  countermeasures  (see  fig.  1). 

The  scenario  defines  an  attack.  The  weapon  sizes,  times  of  burst,  and 
locations  are  specified  for  the  area  being  studied.  The  Locate  Submodel 
produces  attack  environment  data  for  each  unit  area.  The  damage  assessment 
submodel  (ADS)  calculates  the  effects  of  the  weapons  on  the  resources  in 
the  Master  Status  File.  A descriptive  printout  is  then  generated.  The 
"dynamic"  mode  of  TELOS  incorporates  civil  defense  countermeasure  through 
LEMOS.  Civil  defense  operations  determine  the  damage  response  functions 
and  the  resulting  changes  are  expressed  by  improved  states  of  resources 
in  the  Master  Status  File.  Reports  are  generated  during  this  cycle  to  describe 
these  changes  and  the  benefits  derived  from  them. 


3.  Executive  Control 

TELOS  is  controlled  by  an  executive  routine  written  for  the  CDC-3600 
called  the  Generalized  Executive  Control  (GENEC) ; it  is  operable  under  both 
Disk  and  Drum  SCOPE The  control  routine  assumes  that  two  tapes  are 
available:  the  first  tape  contains  all  the  individual  programs  and  the 

control  records  for  the  SORT  utility  that  comprise  TELOS,  while  the  second 
tape  contains  the  parameters  that  may  be  needed  by  the  various  programs 
on  the  first  tape  (see  fig.  2''.  The  use  of  GENEC  is  initiated  by  control 
card  input,  which  tells  GENEC  the  number  of  times  that  the  job  stream  is 
to  be  repeated,  the  number  of  jobs  in  the  stream,  the  position  of  the  first 
job  on  the  program  tape  and  the  position  of  the  first  set  of  parameters  on 
the  parameters  tape.  In  addition,  GENEC  will  be  able  to  pass  five  or  more 
words  of  data  from  one  routine  to  the  next  routine  executed. 

4.  Running  the  Local  Emergency  Operating  System 

The  first  step  in  using  TELOS  with  LEMOS  is  the  selection  of  the 
zone(s)  to  be  studied  and  the  definition  of  the  scenario  to  be  used  in  the 
study  (see  fig.  2). 

In  step  two,  since  TELOS  will  be  controlled  by  GENEC,  the  programs 
comprising  TELOS  (including  those  in  LEMOS)  should  be  stored  on  the  pro- 
gram tape,  and  the  controls  required  by  the  system  should  be  stored  on 
the  parameters  tape.  The  programs  in  LEMOS  require  two  different  types 
of  parameters:  "common"  parameters  and  "specific"  parameters.  The  former 
are  used  by  all  the  programs  in  LEMOS,  and  their  values  are  generally  dependent 
upon  the  scenario  and  the  purpose  of  the  study.  A list  of  these  parameters 
(or  variables)  is  given  in  Table  I.  The  "specific"  parameters  are  used 
to  control  the  functions  of  each  of  the  programs  in  LEMOS.  There  are  three 
"specific"  parameters  which  are  input  to  each  program,  called  RUN-SWITCH, 
TST-SWITCH,  and  PRNT-SWITCH.  The  first  of  these,  RUN-SWITCH,  is  generally 
used  to  select  which  major  functions  a program  should  perform.  The  second 
parameter,  TST-SWITCH,  is  generally  used  for  aid  in  debugging  or  testing  a 
program.  The  last  parameter,  PRNT-SWITCH,  is  generally  used  to  select 
printout  options  within  a program,  e.g.,  which  file(s)  to  list. 

— software  system  operating  on  the  CDC-3600. 


A special  program,  called  the  CONTROL  submodel,  is  being  prepared  to 
run  on  the  CDC-3600.  This  program  will  perform  two  functions:  (i)  accept 

the  original  values  of  both  "common"  and  "specific"  parameters  and  build  the 
relevant  portions  of  a LEMOS  initial  parameters  tape  and  (ii)  allow  changes 
to  the  values  of  the  parameters  during  the  execution  of  the  TELOS  system. 

The  CONTROL  program  will  accept  the  values  of  the  parameters  from  the  parameter 
tape  (batch  mode)  or  from  the  computer  console  (interactive  mode) . 

In  step  three,  a zone  or  local  area  description  is  prepared  in  a 
specific  format.  The  Master  Status  File,  which  contains  the  area's  resources 
except  for  network  data,  is  prepared  on  tape.  A second  tape  is  prepared  con- 
taining network  data. 

In  the  fourth  step,  a job  card  is  prepared  together  with  the  control 
card  described  in  3 above. 

Finally,  in  the  fifth  step  the  cards  and  tapes  are  submitted  for  a 
run  on  the  CDC-3600. 

If  a run  is  to  be  interrupted  at  selected  points  in  the  scenario, 
the  control  cards  reflect  this  decision  and  two  or  more  runs  are  submitted. 


TABLE  I: 


"COMMON"  PARAMETERS  IN  LEMOS 


1.  Hours  from  beginning  of  scenario. 

2.  Duration  of  current  period,  in  hours. 

3.  Sequence  number  of  current  period. 

A.  Flag  to  indicate  whether  or  not  there  was  a previous  period. 

5.  Code  for  minimum  PF  level  for  shelter  spaces  to  be  used. 

6.  Code  for  maximum  height  for  shelter  spaces  to  be  used. 

7.  Low  radiation  level  (RADS)  for  defining  the  basic  operating 
situation  (BOS). 

8.  High  radiation  level  (RADS)  for  defining  BOS. 

9.  Low  fire  level  (fraction  of  area  aflame)  for  defining  BOS. 

10.  High  fire  level  (fraction  of  area  aflame)  for  defining  BOS. 

11.  Codes  for  defining  five  depths  of  debris. 

12.  Zone  number  of  area  being  studied. 

13.  Level  of  PF  provided  by  being  in  automobile. 

14.  Length  of  work  shift,  in  hours. 

15.  Fraction  of  casualties  who  are  ambulatory,  for  fifteen  (15)  injury 
categories . 

16.  Maximum  number  of  teams  to  be  assigned  to  one  problem. 

17.  Identification  of  sanctuary  area  for  this  zone. 

18.  Priority  ranking  of  operations. 

19.  Code  to  indicate  whether  or  not  CD  can  appropriate  resources  from 
residences . 

20.  Code  to  indicate  whether  or  not  population  is  warned. 

21.  Weights  used  to  compute  measure  of  effectiveness. 


A-4 


I 

I 

I 

I 

I 

I 


Appendix  B 

TRANSPORTATION  PROGRAM  DOCUMENTATION 


SECUAiTv  Classification  OF  Tmi*  AaGE  *■*•«  Data  Entaryd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
J BEFORE  COMPLETING  FORM 


«C  T NUM8E4 


44U-873 


2 GOVT  ACC  £ SSION  NO  l ) "fClAlEHTi  CA'ACOG  NUNBtN 


4 T|TLc  (and  Snbntia) 


£)* 


s type  of  REPour  » aeaiqo  COVEAEO 


LOCAL  EMERGENCY  OPERATING  SYSTEM  - LEMOS-  Final  Report  July  1976 


7 . Jo+T**+m*a* 

J.  W./Dunn, 
R.  N. /Hendry 
R. 


--4-  9£A/OJIMinG  QUO  9EPO»t  lUM0Ea 

L-44U-873  _ 

G ° ‘it'  nvjmBC  V 


(T./Lydav 

>o»i»»«;qwcw 


5AHC20-73-C-0253 


V..f  £a#O^M«C  T|QN  same  anO  ACOAESS 

R*esearch  Triangle  Institute 
Post  Office  Box  12194 

Research  Tri ancle  Park,  North  Carolina  27709 


10  pqoGAAm  E.EnEnT  s>«OjF~t  -as* 
A A E A ft  *Os<  JNI  r S i j M a t AS 


DCPA  Work  Unit  41261 




) | 


n Controlling  office  name  a*o  aoo*ES3 

Defense  Civil  Preparedness  Agency 
Washington,  D.C.  20301  


*J.  amh»«T«'9t  A »C«S 

100 


R.  PI 

y-4+f s fpcrt)  If 


>4  MONI  TQAInG  AGEnCv  NAME  ft  AQOAESS fil  dilfarant  irom  Controlling  QUiff 


«S.  SEC  JAITY  C 


. A^/*fS»i  t ' 


Same  as  11. 


Unclassi fied 


•S«.  OEClasSificat'On  OOWnCAaOinG 
schedule 


NA 


'«  OlSTBteu  TION  STATEMENT  (pi  -hit  Raport) 


Approved  for  Public  Release;  Distribution  Unlimited 


17  OiSTAlfuTlON  STATEMENT  t Ol  tha  abtiract  antarad  in  Block  20.  il  dilfarant  from  Raport) 


Same  as  16. 


'•  SuPAL  EMEN  Ti»y  NOTES 


Continuation  of  Earlier  Studies 


t klr  *pAOS  ' Contmua  on  raaataa  nda  if  naeaat 

Civil  Defense 
Local  Operations 
Model ing 

■xecutive  Control 


ary  and  idanttfy  b*  block  numbar) 


Exe< 

Min 

a bVi  ■ 


20  A8\JAACT  (Contlnua  on  ravaraa  aida  If  nac  aaaary  and  idannty  by  block  numbar) 


he  Local  Emergency  Operating  System  (LEM0S)  is  a mul ti -program 
|model  within  a system  of  models  which  are  designed  as  a part  of  a 
'computer  based  system  to  evaluate  local  operating  systems.  The  Defense 
Civil  Preparedness  Agency  Computation  Center  is  a co-developer  with  RTI 
n the  development  of  ADS/LEM0S.  Effort  by  RTI  during  the  past  year 
has  centered  around  the  development  of  the  control  and  transportation  — f 


DD  i 'a""?!  1473  foiriOA  a»  ' «ov  « ii  o«jot.(Ti 


Unclassified. 


SECURITY  CLASSIFICATION  of  This  aaGE  <»***  nata  Fntrraj ) 


SCCU"ITV  Ct  AJS**'C  *t'On  Of  Twi*  0*<«  Bnffd) 


submodels  and  at  the  same  time  has  continued  to  improve  the  procedures 
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This  report  describes  the  essential  features  of  the  control  and 
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to  development  but  with  the  specific  objective  of  gaining  support  for 
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Adequate  developmental  testing  is  recommended  before  any 
demonstrations  of  the  planning  role  are  undertaken. 
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TRANSPORTATION  SUBMODEL 

Section  I:  Abstract  and  Run  Description 

Abstract 

Inputs  to  the  Transportation  Submodel  comprise  the  Links  File  and 
Control  File,  both  prepared  manually  by  the  user,  and  the  Resource  File  which 
is  generated  by  the  Problem  Definition  Submodel.  The  Links  File  describes 
the  network  in  each  Unit  Area  in  its  initial  state.  The  Control  File  is  a 
list  of  Unit  Areas,  possibly  a subset  of  the  total,  selected  for  processing. 

The  Resource  File  indicates  changes  in  state  of  each  Unit  Area  via  the  data 
variable  called  Basic  Operating  Status,  BOS. 

The  program  computes  the  shortest  paths  between  all  nodes  in  a Unit 
Area  network  and  between  all  Unit  Areas. 

The  primary  output  is  the  TVL-REC  file,  a list  of  the  minimum  travel 
times  between  all  Unit  Areas  defined  in  the  resource  file.  Other  outputs 
are:  print-outs  of  the  distance  and  policy  matrices  before  and  after  the 
shortest  paths  are  computed;  various  working  or  temporary  files;  and  files, 
either  temporary  or  permanent,  of  the  distance  and  policy  matrices,  accor- 
ding to  the  user's  desires  and  available  computer  installation  options. 

Run  Description 

History 

In  order  to  affect  civil  defense  countermeasures,  resources,  equipment, 
and  personnel  must  be  dispatched  to  areas  requiring  remedial  action  in  time 
to  be  of  benefit.  The  need  for  knowing  the  route  of  minimum  travel  time 
between  the  origin  of  the  resources  and  the  destination  in  need  of  remedial 
action  is  self-apparent.  Also  apparent  is  the  need  to  choose  between  competing 
demands  and  alternate  available  resources.  The  solution  of  this  complex 
allocation  problem  requires  knowledge  of  the  relative  minimum  travel  times  to 
dispatch  resources  to  demands. 

Specifications  and  Memoranda 

The  program  is  written  in  FORTRAN  IV.  It  has  been  compiled  in  FORTRAN-H, 
link  edited,  and  executed  on  an  IBM  370/175  computer  using  0S/360,  Level  21.6. 
While  every  attempt  was  made  to  avoid  IBM-specific  features  of  the  language 
in  order  to  permit  use  of  the  program  on  other  machines  with  FORTRAN  compilers. 


: 


some  changes  may  still  be  necessary  especially  in  regard  to  I/O  statements. 

An  estimate  of  available  computer  core  memory  at  the  ultimate  user's 
installation  indicated  that  the  progra.i  be  written  so  that  in  execution  it 
would  require  not  more  than  250  x 210  bytes  (1  byte  = 8 bits)  on  the  IBM  370. 
This  specification  has  been  met  exclusive  of  I/O  buffers  and  out-of-core 
(disk  or  tape)  data  storage  requirements. 

Procedure 

See  Section  5,  Operating  Instructions 

Method 

The  program  uses  Floyd's  method  for  determining  the  paths  of  shortest 
distance  between  all  nodes  in  a network.  See  "An  Appraisal  of  Some  Shortest- 
Path  Algorithms"  by  Stuart  E.  Dreyfus  in  ORSA,  Vol . 17  #3,  1969. 

Formulae  Used 

Given  the  n x n initial  matrix  of  directly  linked  nodes  in  a network: 

D»  = (d?J) 
d?,i 

d?  . = °°  if  nodes  i and  j are  not 
directly  linked, 

d?  . < °°  if  nodes  i and  j are 
directly  linked. 


where 


and 


Floyd's  algorithm  computes  successively  for  k = 0,1,2, ,n 


m1n  <d?,r  + dk,o> 


i . j = 1,2, ,n . 


Restrictions 

There  are  restrictions  imposed  by  FORTRAN  dimension  statements: 

75  = max.  number  of  nodes  in  unit  area  network 
400  = max.  number  of  unit  areas. 

These  limitations  may  be  changed  by  changing  the  dimension  statements. 
Options 

Out  of  the  total  number  of  unit  area  networks  represented  in  the  Links 
File,  only  those  specified  in  the  Control  File  will  be  processed. 

Device  numbers  for  files  read  from  IN  and  written  on  I0UT  may  be  changed 
in  the  BLOCK  DATA  subprogram. 


The  Control  File  is  read  from  device  IN.  Print-out  reports  are  written 
to  device  IOUT.  Other  data  files  use  device  numbers  set  in  the  program.  A 
complete  description  is  in  Section  3,  Input-Output  Description. 

Accuracy 

The  distances  between  nodes  are  rounded  to  the  nearest  tenth  of  a mile. 
Computed  averages  of  these  distances  are  rounded  to  the  nearest  tenth  of  a 
mile. 
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Anticipated  Usage 

The  Transportation  Submodel  is  one  within  the  mul ti-program  model  Local 
Emergency  Operating  System  (LEMOS),  which  is  to  be  integrated  into  the  Test 
and  Evaluation  of  Local  Operating  Systems  (TELOS)  model,  which,  in  turn,  is 
to  function  under  the  control  of  an  executive  routine,  GENEC. 

The  Transportation  Submodel  comprises  two  major  steps.  The  first  prepares 
a data  file  of  distances  between  unit  areas  for  use  by  the  second.  It  is 
anticipated  that  the  first  will  be  used  as  a stand-alone  program,  even  possi- 
bly outside  GENEC.  The  transportation  network  data  base  (Links  File)  is 
developed  by  Unit  Areas  over  a period  of  time  for  a geographical  study  area. 
Each  Unit  Area  network  as  it  is  developed  can  be  processed  by  the  first 
program  and  the  results  analyzed  off  line  for  consistency  by  comparison  with 
the  data  source  documents,  topographic  or  highway  maps.  Subject  to  any 
changes  or  corrections,  the  verified  results  may  then  be  added  to  the  data 
file  input  for  the  second  program.  The  second  step  computes  shortest  paths 
between  unit  areas.  It  is  anticipated  that  it  will  be  used  in  a time-itera- 
tive scenario  under  the  control  of  GENEC. 

Timing  Factors 

Execution  time  in  seconds  by  CPU  and  I/O  for  the  case  of  8 unit  areas 
containing  a total  of  120  nodes  and  from  6 to  26  nodes  per  unit  area  was: 


Actual  Time  in  Seconds 


CPU 

1/0 

1st  Step 

(find  all  shortest 

paths  in  all  unit  areas) 
2nd  Step 

(find  all  shortest 

10.2 

28.5 

oaths  between  all  unit 
areas ) 

0.6 

9.5 

To  obtain  timing  estimates  of  a larger  scale  problem,  say  400  unit  areas, 
the  above  times  multiplied  by  50  (=  400/8)  are  interpreted  subject  to  the 
following  considerations. 

First  of  all  the  technique  of  simple  one-point  extrapolation  is  question- 
able, but  it  is  the  only  data  available. 

The  1/0  times  include  time  required  to  buffer  the  print-out  of  all 
distance  and  policy  matrices  before  and  after  execution  of  the  shortest  paths 
algorithm.  Actual  printer  time  is  not  included.  These  print-outs  would 
not  occur  in  regular  usage. 

The  CPU  and  1/0  times  above  were  obtained  on  an  IBM  370  operated  under 
MVT  (multiprogramming  with  a variable  number  of  tasks.)  The  number  of  other 
jobs  and  the  job  mix  in  the  system  has  some  effect  on  these  times.  In  a 
different  environment  or  on  a different  computer  dedicated  to  the  task,  differ- 
ent results  would  occur. 

The  following  time  estimates  for  a 400  unit  area  problem  are  felt  to  be 
overly  pessimistic,  except  for  the  second  step  1/0  time  which  is  probably  too 
low. 

Estimated  Time  in  Seconds 


CPU 

I/O 

1st  Step 

510 

1,425 

2nd  Step 

30 

475 

Section  II:  Flow  Charts 


Narrative  descriptions  in  lieu  of  flow  diagrams  are  provided  for  the 
following  utility  type  programs. 

1.  Subroutine  RANK 

Subroutine  RANK  is  used  to  rank  order  the  elements  of  a vector 
array  in  ascending  order.  In  particular.it  is  used  in  NETWRK  to  rank  the 
node  identification  numbers  of  a given  unit  area  network.  The  purposes  of 
rank-ordering  are  to  facilitate  inspection  of  the  results  of  data  preparation 
and  shortest  paths  computations  and  to  determine  the  method  of  finding  equal 
node  i.d.'s  in  two  unit  area  networks.  The  vector  of  unranked  values  is 
input  to  the  subroutine;  the  ranked  values  are  returned  in  a second  vector. 

A third  vector  associates  the  new  index  of  a ranked  value  with  its  original 
index  before  ranking.  Briefly,  to  illustrate,  consider  unranked  values 
(c,b,a)  returned  as  (a,b,c): 


Original  Order  Value 

1 c 

2 b 

3 a 


Ranked  Order  Value 

3 a 

2 b 

1 c 


Thus, the  elements  of  the  third  array,  1 1, are: 


11(1)  = 3,  11(2)  = 2,  11(3)  = 1. 


2. 

ei ther 


Functions  COLSUM  and  ROWSUM 

These  two  functions  operate  on  a given  matrix,  (d.  .),to  compute 

' * J 


or 


COLSUM  = 
ROWSUM  = 


n 


d 


i J 


In  particular,  these  functions  are  used  in  ALLNET  to  compute  the  total 
distance  from  all  origins,  i , to  a specified  destination,  j,  or  from  a speci- 
fied origin,  i,  to  all  destinations,  j,  where  (d.  .)  is  the  matrix  of  shortest 
distances  between  all  i and  j.  In  either  case  the  specified  node  is  a boun- 
dary node  between  two  unit  areas. 


3.  Function  STORE  and  Function  UNMASK. 


Function  STORE  and  its  entry  type  function,  UNMASK,  are  used 
respectively  to  compress  and  extract  the  five  hierarchial  identifications  of 
a unit  area  network  in  a single  word  as  summarized  in  the  following  table. 


Sing' 

1 e-word  bit  configuration 

Magni tude 

Identification  Vector  1 

1 

2 

3 

4 

5 

Zone 

HIARCH ( 1 ) 

///// 

25-1 

EOC 

H I ARCH ( 2 ) 

mm 

26-l 

Group 

H I ARCH ( 3 ) 

mm 

26-l 

Sector 

HIARCH(4) 

miii 

26-l 

Unit  Area 

HIARCH(5) 

iiiiinn 

29-l 

Total  bits 

5 

6 

6 

6 

9 

32 

Incl usi ve 

Positions 

32-28 

27-22 

21-16 

25-10 

9-1 

To  illustrate,  the  EOC  identification  of  a unit  area  stored  in  the  array 
HIARCH(2),  is  stored  in  bit  positions  27-22  in  the  single  word  NETS  (I)  for  a 
specified  Ith  network.  This  is  accomplished  by  multiplying  HIARCH  (2)  by  221 
(which  moves  the  value  to  the  proper  bit-field)  and  then  computing  the  logical 
sum  of  the  result  and  NETS  (I).  To  extract  the  EOC  from  NETS  (I),  the  logical 
product  of  NETS  (I)  and  the  2nd.  bit-field  configuration  above  is  computed  and 
then  the  result,  in  bit  positions  27-22,  is  divided  arithmetically  by  221  to 
move  it  to  the  lower  order  bit  positions,  1-6. 

The  masks  (bit-field  configurations)  and  powers  of  2 are  computed  by  a 
first  call  to  STORE.  For  machines  with  less  than  32-bit  words,  these  values 
must  be  re-coded.  For  machines  with  greater  than  32-bit  words,  re-coding  is 
unnecessary  but  it  may  be  desirable  in  order  to  increase  the  maximum  permis- 
sible magnitudes. 

4.  Subroutine  MATCH 

This  subroutine  compares  arithmetically  the  elements  of  two  integer 
arrays  for  equality  between  any  elements.  It  returns  a code  of  either  0 or  1 
indicating  either  "no"  or  "yes".  If  there  is  equality,  the  indices  of  the 
equal  elements  are  returned.  The  subroutine  is  used  in  subroutine  SETUP  to 
create  the  vector  of  node  identification  numbers  from  the  list  of  nodes  con- 
tained in  the  Links  File. 


SUBROUTINE  AND  FILE  POSITIONS  IN  TRANSPORTATION  SUBMODEL 


Create  matrix  of  distances 
between  directly  linked 
nodes  in  each  unit  area 


NETWORK 
*|  Mai  n ] 

—j Block  Data 
— jCntlrd 
-H  Setup 


|Fi 1 e In,  Control  File  j 
fFTTe"  51,  Links  File  I 


Compute  shortest  paths 
between  nodes  in  unit 
area 

File  30,  # networks, 
packed  network  i .d. 
numbers 


Rank  ■ j 
^SpathJ 


Prnout 


Unmask 


Distance 
and  Policy 
Matrices 


File  53;  Network  i.d.  #, 
# nodes,  node  id's, 
matrix  of  shortest  paths 


Compute  distances  between 
directly  linked  unit  areas 


ALLNET 

Main  r 

■j  Rows  urn 
■f  Col  sum 


Netprn 


Working  Files: 
/'File  70 
File  71 
File  80 
^File  81 
Data  Files: 


File  91 ; distance  from 
A to  boundary 


File  92;  distance  from 
boundary  to  B 


Time  Iteration  Scenario 


t.  J*. 


SUBROUTINE  AND  FILE  POSITIONS  IN  TRANSPORTATION  SUBMOOEL  (Continued) 


Damage  assessment 


Resource  File 


Compute  shortest  paths 
between  unit  areas 


Format  the  data  for 
TVL-REC  File 


Y 


J 


NETWORK  2 of  6 


Section  III:  Input-Output  Description 

The  general  order  of  arrangement,  creation,  and  utilization  of  the  18 
files  in  the  Transportation  Submodel  is  shown  in  Section  II,  Flow  Diagrams. 
Details  concerning  the  preparation  and  format  of  the  two  manually-prepared 
files.  Control  and  Links,  are  contained  in  the  main  report.  This  section 
will  primarily  describe  the  requirements  of  the  submodel -generated  files. 

All  files  are  sequential  files.  In  the  testing  and  evaluation  of  the 
Transportation  Submodel,  all  files  except  IOUT  (the  printer)  were  maintained 
on  on-line  disk  packs.  The  REWIND  instruction  occurs  numerous  times  for 
closing  files  not  in  use  to  conserve  computer  core  (used  for  I/O  buffers) 
as  well  as  to  reposition  files  for  repetitive  processing.  The  file,  IN,  which 
is  normally  the  system  input  device,  i.e.  card  reader,  is  also  rewound;  this 
could  be  prohibitive  at  certain  computer  installations  if  the  device  is  not 
indeed  a disk  or  tape. 

Generally,  all  files  are  unformatted,  exceptions  being  the  files  used 
for  prepared  input  (IN,  the  Links  File,  and  the  Resource  File)  and  formal 
output  (on  device  IOUT  and  the  TVL-REC  File.)  Samples  of  the  formatted  data 
files  are  shown  in  Section  IV,  Test  Data. 


Logical  No.  records  as  function 


Name 

Device 

Format 

Record  Length  Blocksize 

of 

no.  unit  areas,  n 

Control 

IN 

FB 

80 

12960 

n/7i/ 

Links 

51 

FB 

700 

13000 

n.a. 

53 

VBS 

13022 

13026 

n 2/ 

IOUT 

IOUT 

FBA 

132 

132 

n.a. 

30 

VBS 

1608 

11260 

1 

70 

VBS 

1608 

11260 

n 

71 

VBS 

1608 

11260 

n 

80 

VBS 

1608 

11260 

n 

81 

VBS 

1608 

11260 

n 

TOBOUN 

91 

VBS 

1608 

11260 

n 

FRMBOU 

92 

VBS 

1608 

11260 

n 3/ 

Resource 

61 

FB 

129 

2580 

n.a. 

93 

VBS 

1608 

11260 

n 

60 

VBS 

1608 

11260 

n 

70 

VBS 

1608 

11260 

n 

80 

VBS 

1608 

11260 

n 

55 

VBS 

1608 

11260 

n 

TVL-REC 

94 

FB 

20 

13020 

n2 

1/ 

1 record  per  link  for  each  unit  area  network 

+ n '999'  separator 

records 

2/ 

Will 

vary  with 

the  type  print-out  being  produced 

3/ 

Depends  on  number  of  Land  Use  Categories 

in 

each  unit  area. 

Suninary  of  Data  Files  Used 
in  Transportation  Submodel 


LINKS  FILE  Record  Format-7 


COBOL 

Variable 

COBOL 

Format 

Fortran  Fortran 

Var iable(s)  Format 

Card 
Col umns 

REMARKS 

ZONE 

9 

H I ARCH ( 1 ) 

11 

1 

Zone  i.d.  number 

EOC 

9 

H I ARCH ( 2 ) 

11 

2 

EOC  i.d.  number 

GROUP 

99 

H I ARCH ( 3 ) 

12 

3-4 

Group  i.d.  number 

SECTOR 

99 

HIARCH(4) 

12 

5-6 

Sector  i.d.  number 

NETWORK-NO-7 

999 

HIARCH(5) 

13 

7-9 

Unit  Area  i.d.  number 

FWD-NODE 

999 

IFORW 

13 

10-12 

Forward  node  i.d.  number 

none 

X 

LVLFOR 

11 

13 

Level  of  forward  node 

BWD-NODE 

999 

I BACK 

13 

14-16 

Backward  node  i.d.  number 

none 

X 

LVL8AC 

11 

17 

Level  of  backward  node 

LINK-NO 

9999 

LINKWC 

14 

18-21 

Link  i.d.  number 

LINK-NAME 

x(ll) 

not  used  , 

22-32 

Link  name 

UA-NO-L 

999 

not  used 

33-35 

Uni t area , left  of  link 

NTWK-NO-L 

999 

not  used 

» 23X 

36-38 

Network,  left  of  link 

UA-NO-R 

999 

not  used] 

39-41 

Unit  area,  right  of  link 

NTWK-NO-R 

999 

not  used  1 

42-44 

Network,  right  of  link 

LENGTH 

9V9 

AMILES 

F3.1 

45-47 

Length  of  link  in  miles 

WIDTH 

9 

not  used 

IX 

48 

Width  of  link 

TYPE 

XX 

TYPE 

A1 

49 

Use- type  of  link 

ITYPE 

11 

50 

0 / 

Direction  of  link 

none 

not  used 

50-100— 

Room  for  further  expan- 

sion;  variables  for 
queuing  model . 


Notes: 

1/  This  file  should  be  sorted  by  ZONE,  EOC,  GROUP,  SECTOR,  and 
NETWORK-NO  in  ascending  order. 

2 / A record  with  NETWORK-NO  - 999  separates  networks. 

2J  For  development,  LRECL  = 80  and  columns  73-80  contained  sequential 
card  identification  numbers. 


Description  of 
Links  File  Records 


Interior  to  unit  area;  not  shared  with  other  areas 

Shared  between  unit  areas 

Shared  between  sectors 

Shared  Between  groups 

Shared  between  EOC's 

On  the  extreme  boundary  of  zones 


Codes  Used  in 
Links  File 


Control  Cards 

Format 

Format 

Col umns 

Variable  name 

13 

1-3 

NETNUM(I);  i.d.  of  unit  area  to 
be  processed 

11 

4 

IFLAG  - no  longer  used 

13 

5-7 

I OR I G - no  longer  used 

13 

8-10 

IDEST  - no  longer  used 

The  above  format  occurs  7 times  per  card. 

A value  of  NETNUM(I)  <_  0 terminates  the  string  of  network  numbers 
This  format  is  used  in  Subroutine  CNTLRD. 


Description  of 
Control  File  Records 


Test 

Data 


Section  IV:  Test  Data 


Test  data  were  derived  from  maps  of  the  Detroit  area  for  8 unit  areas. 
Roads  of  interest  were  manually  transcribed  from  published  maps  to  a working 
sheet.  Link  names  and  identification  numbers  for  nodes,  links,  and  unit 
areas  were  written  on  the  working  sheet.  These  data  were  transcribed  to 
coding  forms  in  accordance  with  the  Links  File  format  described  in  the  main 
report  with  the  exception  of  the  inclusion  of  sequential  card  numbers  in 
columns  73-80.  Listings  of  the  Links  File  and  the  Control  File  follow. 

For  purposes  of  development  and  test,  an  existing  Resource  File  was 
modified  by  truncating  the  records  to  80  columns.  To  this  modified  file  were 
added  records  for  LUC  = 5 and  for  other  unit  areas  in  the  transportation 
network  but  not  in  the  Resource  File.  A listing  of  the  modified  Resource 
File  follows. 

Printed  results  before  and  after  the  shortest  paths  computation  using 
the  Control,  Links,  and  Resource  data  files  as  input  follow.  Comparison  of 
the  distance  and  policy  matrices  with  the  network  as  drawn  on  the  working 
sheet  is  a manual/visual  operation.  Selecting  key  links,  nodes  at  extreme 
f points  and  measuring  some  distances  on  the  map  helped  uncover  errors  in  the 

Links  File,  aided  in  de-bugging  the  program  and  verified  the  accuracy  of  the 
final  results. 

There  were  few  predetermined  results  that  could  be  used  except  for 
shortest  paths  using  oodward  Avenue  or  the  Chrysler  Freeway. 


>• 
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400014090191000000000000007000007000007000001000 OOOOOOOOOOO p 00000000000000000280 
4000  141  00  19 100000y0o000000700000700000 700000 1000000  0 0000000000000000000000000270 
4000141  101  91  001  390  0 73 360263670 0070002e 36 7 0 0 0 1 0 0 0 0 0 l 0 00 0 v 00 0 0 00 0 0 0 0 0 0000 0 0 00 0 0280 
4000l4l2019lo000o0000000007000007000007000y01000000c00y0y000000000000y0000c00290 
400014130l910o0000000u0000400000400000400000100c0000000000u000u0000000000000c300 
40  0 0 141  70  19  l OoOOO 0000000005000005000005000001 0000O0 00000 00U0P0OO00000 0000000031  0 
40001422019100071 00000y0o04000004000y0400000100y0yl00000000000000000000000000 320 
4000 142401 9i Oooooooyoooooosooooo 30 0000300000 ioooopooooooooooooooooooooooyoono33y 
4000l4330191000000000000002000y02000002000001o0u000000o0000000000000j00000u00540 
4000 14390 19 100000o00000000200000200000200000100o00000000000y00y0000000000j000380 
4000  144 00  191  000000 00 0000002O 0000 20 000020y000ly0  0000o000000000000y000000 0 00000380 
90001441019100000000000000700000700000700000100000000000000000000000 000000000370 
4OO01446O141o0O0O0000Oy000700O007O0OyO7000OOlOOo00O0CO0000OO0O0OO000000O0CUO039G 
400014480191 00040 0036650 12575000 700 01 25 750001000001 00000000 00 0000000000000000390 
40001449019  10 8 004000466001 62 30 0070000  16230001  0000010000000000000000000000000040o 
40001 452019 100000000000000700000700000700000100000000000000000000000000000000 41 U 
4000l453019loo0000000y00007y00y0700y007oy00010000000000000000000000y00000y000420 
400014590 191 000 12001 0000007000007u00u0700000l 0000751 0000000000000000000o000o0430 
4000i4oiui9iy002d0000i000i2000yi20000i2y0000iyycy0i0000000y0000000y0000y00o00440 
4oyoi463oi9ioy0yooooooyyo040yooo4oyoyy40oooyooooooccooyyoyocooooyouooyoooooyo45o 


2105 

2905 

3105 

7605 

8905 


1 
1 
1 
1 
l 
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00000460 

00000470 

00000480 

00000490 

00000500 
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. "W. 


1 112 

3 

3 

153 

142 

7LlvERNOIS 

21 

21 

3 

31.0 

SI 

OOOOOO  1 0 

1112 

3 

3 

383 

153 

7t  IGm  T MI  RO 

76 

7e 

3 

30.8 

SI 

00000020 

1112 

3 

3 

763 

363 

6E  I Gn  T MI  SO 

76 

76 

3 

30 .6 

SI 

00000030 

1112 

3 

3 

uu  3 

7d  3 

5EIGHT  MI  SO 

89 

89 

3 

30.8 

SI 

00000040 

1 112 

3 

3 

555 

443 

4E IGm T MI  SO 

89 

89 

3 

3o . e 

SI 

00000050 

1 112 

3 

3 

624 

555 

oCQNA.i 7 4 VE 

29 

29 

3 

31 . 1 

SI 

00000060 

1112 

3 

3 

624 

544 

4S£V£m  mi  so 

3 

3 

29 

290.5 

SI 

00000070 

1112 

3 

3 

544 

534 

40EUUINUME 

3 

3 

29 

290.7 

SI 

00000080 

1112 

3 

3 

534 

524 

3DE Du I NOSE 

3 

3 

14 

140.3 

SI 

00000090 

1112 

3 

3 

524 

504 

OMCMICHULS  9 

3 

3 

14 

140.2 

SI 

00000100 

1 112 

3 

3 

504 

414 

7MCNICHOLS  * 

3 

3 

4 

40,4 

SI 

00000110 

1112 

3 

3 

4 14 

3o4 

8MCNICHULS  4 

3 

3 

4 

40,6 

31 

00000120 

1112 

3 

3 

364 

334 

9MCNICMULS  P 

3 

3 

4 

40.5 

31 

00000130 

1112 

3 

3 

334 

834 

IOMCNICmOLS  S 

3 

3 

4 

40.3 

SI 

00000140 

1112 

3 

3 

634 

142 

IPASKSIDE 

3 

3 

21 

211.9 

FI 

00000150 

1112 

3 

3 

371 

142 

7S£v£\  mi  RQ 

3 

3 

3 

31.3 

St 

00000160 

1112 

3 

3 

431 

371 

68EVEN  MI  RO 

3 

3 

3 

31.0 

SI 

00000170 

1 112 

3 

3 

544 

431 

5SEVE,'»  mi  so 

3 

3 

3 

30.6 

SI 

00000180 

1112 

3 

3 

371 

363 

umJODAAPD  A v 

3 

3 

3 

31.1 

SI 

00000190 

1112 

3 

3 

364 

371 

SmJOOnARC  Av 

3 

3 

3 

31.1 

SI 

00000200 

1112 

3 

3 

443 

431 

5CmRY5LEP  fy 

3 

3 

3 

31.0 

S4 

00000210 

1112 

3 

3 

431 

421 

uCHKYSLEP  fr 

3 

3 

3 

30.2 

S4 

00000220 

1112 

5 i 

504 

421 

SCmpyslER  Fr 

3 

3 

3 

30,9 

Su 

00000230 

00000240 

1123 

a 

4 

414 

364 

6MCNICH0LS  « 

3 

3 

4 

40.6 

SI 

00000260 

1123 

a 

4 

364 

334 

9MCMICH0LS  4 

3 

3 

4 

40.5 

SI 

00000270 

1123 

a 

4 

334 

834 

lOMCMCHOLS  9 

3 

3 

4 

40.3 

SI 

00000280 

1 123 

a 

4 

292 

634 

1THUMS0N 

4 

4 

31 

311.5 

SI 

00000290 

1123 

a 

4 

285 

292 

2J. LODGE  For 

4 

4 

31 

310.3 

54 

00000300 

1 123 

a 

4 

275 

285 

1 J •LODGE  Fay 

4 

4 

0.4 

S 4 

00000310 

1123 

a 

4 

315 

275 

37 UxEOO 

4 

4 

0.2 

SI 

00000320 

1123 

4 

4 

3u5 

315 

2TUXEQ0 

4 

4 

0.5 

SI 

00000330 

1123 

a 

4 

395 

345 

1 7UXED0 

4 

4 

0.5 

SI 

00000340 

1123 

4 

4 

321 

398 

8UAVISUN  EXP 

4 

4 

4 

40.2 

Su 

00000350 

1123 

4 

4 

351 

321 

7DAv ISOM  ExP 

4 

4 

4 

40.5 

S 4 

00000360 

1123 

4 

4 

401 

351 

60  A V I SOi<  EXP 

4 

4 

4 

40,5 

Su 

00000370 

1123 

4 

4 

663 

401 

50AVISUIY  EXP 

4 

4 

4 

40.3 

Su 

00000380 

1123 

4 

4 

315 

321 

IHAMILTOy  Av 

4 

4 

4 

40.7 

SI 

00000390 

1 123 

4 

4 

321 

334 

2MAMILTON  AV 

4 

4 

4 

41.4 

SI 

00000600 

1123 

4 

4 

345 

351 

lnUODnASO  av 

4 

4 

4 

40.7 

SI 

00000410 

1123 

4 

4 

351 

364 

2a JGOnAPO  Av 

4 

4 

4 

41.1 

31 

00000420 

1 123 

4 

4 

395 

401 

l OAKLAND  Avt 

4 

4 

4 

40.7 

SI 

0 0000430 

1123 

4 4 
999 

40  1 

4 1 4 

20AKI.A.VD  AvE 

4 

4 

4 

40.9 

SI 

00000440 

00000450 

1 1 2<a  i a 

14 

524 

50  4 

6MCMCH0LS  R 

3 

3 

14 

140.2 

SI 

00000460 

1 12«*t 

4 

14 

534 

524 

SOEGUIaOHE 

3 

3 

1« 

140.3 

SI 

00000470 

1 12**  1** 

14 

612 

534 

1 M IMNt SO  T A 

29 

29 

14 

140.3 

SI 

00000400 

1 12«lu 

14 

632 

612 

IE.  9EVA0A 

29 

29 

14 

140.5 

SI 

OOOOOuPO 

1 12uiu 

14 

692 

oS2 

2E.  4EVA0A 

29 

29 

14 

1«1  .0 

SI 

00000500 

1124  14 

14 

735 

692 

3E.  NEVAOA 

29 

29 

14 

141.0 

SI 

00000510 
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1 1 2 a i a 

14 

715 

725 

IVAN  DYKE 

14 

14 

0.5 

SI 

00000520 

t 12*414 

14 

725 

815 

lwCMCHULS  t» 

14 

14 

0.8 

SI 

00000530 

112414 

14 

915 

805 

mt  elliut r 

14 

14 

0.9 

SI 

00000540 

112414 

1« 

095 

975 

1C AN  IFF  Avfc 

14 

14 

0.3 

SI 

00000550 

112414 

14 

o 75 

875 

2 C A N I F F AVE 

14 

14 

0.3 

SI 

OOOOOSoO 

112414 

14 

9 75 

575 

1CANIFF  avc 

14 

14 

u.3 

51 

00000570 

112414 

14 

575 

585 

ICQNAiOr  AVE 

14 

14 

0.5 

SI 

00000580 

112414 

14 

585 

485 

1CA8PENTES 

14 

14 

0.8 

SI 

00000590 

112414 

14 

485 

495 

2CAAPE'9rt4 

14 

14 

0.3 

SI 

00000800 

112414 

14 

485 

511 

1 ChP  Y SLEP  F 1 

14 

14 

14 

140.5 

S 4 

00000910 

112414 

14 

511 

9oi 

30AVISUH  Exp 

tu 

14 

14 

140.  i 

S 4 

00000920 

112414 

14 

511 

504 

2C8MYSLES  FY 

14 

14 

14 

140.8 

S 4 

00000630 

112414 

14 

eOl 

524 

5MCN I CHOLS  p 

14 

14 

14 

141.0 

SI 

00000940 

112414 

14 

921 

oO  1 

4MC*ICrlOLS  4 

14 

14 

14 

140.5 

SI 

00000850 

112414 

14 

981 

821 

JMCMICHGLS  4 

14 

14 

14 

140.5 

SI 

00000860 

112414 

14 

915 

981 

24CMICH0LS  9 

14 

14 

14 

140.3 

SI 

00000970 

112414 

14 

591 

511 

2'JAvISO'y  cap 

14 

14 

14 

iuo.8 

S4 

00000680 

112414 

14 

821 

591 

ID AV ISON  E. 

1« 

14 

14 

140.5 

SI 

00000990 

112414 

14 

90  1 

o 1 2 

4car.A,»r  ave 

14 

14 

t« 

140.3 

SI 

00000  700 

112414 

14 

591 

oO  1 

3ClmAnt  AVt 

14 

14 

14 

140.2 

SI 

0000071  0 

112414 

14 

585 

591 

2CUNA,yI  ave 

14 

14 

l« 

140.5 

SI 

00000720 

112414 

14 

ol2 

90  l 

19YA.N  STPEE  7 

14 

14 

14 

140.5 

SI 

00000710 

112414 

14 

o92 

691 

1 •tODNO  KUAO 

14 

14 

14 

140.5 

sv 

00000740 

112414 

14 

o81 

o 75 

2*40040  PUAD 

14 

14 

14 

141.0 

SI 

00000750 

999 

00000790 

111221 

21 

905 

105 

4F IkELANE 

21 

211.0 

FI 

00000770 

111221 

21 

795 

905 

5F ISEuAnE 

21 

211.0 

FI 

00000790 

1 1 1221 

21 

75 

795 

OEI&mT  «I  90 

21 

210.5 

SI 

00000790 

111221 

21 

151 

75 

SEllirtT  <1  PU 

79 

79 

21 

211.0 

SI 

00000800 

111221 

21 

151 

142 

7l  I yEPNO IS 

21 

21 

3 

31.0 

31 

00000310 

111221 

21 

814 

1h2 

lPAHKSIOE 

1 

3 

21 

211.9 

FI 

00000820 

1 1 1221 

21 

314 

114 

1 l">Cfi  IChOLS  4 

21 

21 

21 

210.9 

SI 

00000830 

1 1 1221 

21 

114 

184 

5tl vEpnQI  S 

21 

21 

31 

310.5 

SI 

00000940 

111221 

21 

124 

44 

I PUP  I 1AM 

21 

21 

31 

311.4 

SI 

00000850 

111221 

21 

44 

105 

6CUUZENS  F.YY 

21 

21 

31 

311.1 

S 4 

00000890 

1 1 1221 

21 

51 

44 

4*V0MI;iG  AyE 

21 

21 

21 

210.9 

SI 

00O00870 

111221 

21 

ol 

51 

5.y YOMI.yG  avE 

21 

21 

21 

211.0 

SI 

00000880 

1 1 1221 

21 

75 

91 

9«Y0"4ING  Av£ 

21 

21 

21 

211.0 

St 

00000890 

111221 

21 

142 

114 

6L I V E 9N0 1 S 

21 

21 

21 

211.0 

SI 

00000900 

111221 

21 

91 

805 

PSElEN  "l  40 

21 

21 

21 

210.5 

SI 

00000910 

111221 

21 

142 

61 

33tVtN  4 I PD 

21 

21 

21 

210.5 

SI 

30000920 

111221 

21 

51 

105 

11MCNICM0LS  R 

21 

21 

21 

211.0 

SI 

00000910 

1 1 1221 

21 

114 

51 

1 2PCN ICHULS  » 

21 

21 

21 

211.0 

SI 

00000940 

949 

00000950 

1 12429 

29 

912 

514 

I^INpESOT  A 

24 

29 

14 

140.8 

St 

000009sO 

112429 

29 

o 12 

612 

It.  NEVADA 

24 

29 

14 

t “0 . 5 

SI 

00000970 

l 12429 

29 

992 

912 

26.  i £ v A 0 A 

24 

29 

14 

141.0 

SI 

00000990 

1 12429 

29 

715 

692 

JE.  NEVADA 

24 

29 

14 

141.0 

31 

00000940 

1 12429 

24 

745 

715 

2v an  OyaE 

24 

24 

0.5 

SI 

00001000 

1 12429 

29 

755 

745 

3VA.4  OYXE 

24 

24 

1.0 

SI 

00001010 
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I 

I 

I 


112424 

24 

755 

715 

lEIGrlT  Ml  40 

24 

241.0 

SI 

00001020 

1 12424 

24 

715 

o55 

2EI&"7  -*1  RQ 

24 

241.0 

SI 

ooaoioio 

1 12424 

24 

o55 

555 

SEIGnT  41  40 

24 

241.0 

SI 

00001040 

1 124^4 

24 

0 2 4 

555 

bCUivant  a vE 

24 

24 

3 

31  . 1 

SI 

00001050 

1 12424 

24 

824 

544 

4SEv£<  41  RO 

3 

3 

24 

240.5 

SI 

00001 OeO 

1 12424 

24 

544 

534 

4y£'JUl.04£ 

3 

3 

24 

240.7 

SI 

00001070 

1 12424 

24 

812 

824 

5CJNA.VT  avE 

24 

24 

24 

240.8 

SI 

00001 080 

1 12424 

24 

841 

824 

3 S E v E N Ml  w to 

24 

24 

24 

2 40.5 

SI 

00001040 

1 12424 

24 

701 

641 

2SEvEn  '4 1 80 

24 

24 

24 

241.0 

SI 

00001 100 

1 12424 

24 

745 

701 

ISEVEi  4 1 fio 

24 

24 

24 

241.0 

SI 

0000  1 1 1 0 

1 12424 

24 

84  1 

632 

2HYA.V  street 

24 

24 

24 

240.5 

SI 

aoooi 120 

1 12424 

24 

855 

84  1 

3HTAI,  STREET 

24 

24 

24 

241.0 

SI 

OOOO 1 1 io 

1 12424 

24 

701 

842 

4MUON0  ROAD 

24 

24 

24 

240.5 

SI 

00001140 

l 12424 

24 

715 

701 

5MOIA0  ROAO 

24 

24 

24 

241.0 

SI 

00001150 

999 

00001180 

112331 

31 

855 

885 

1 ME Y £ W S ROAD 

31 

310.5 

SI 

0000 l 170 

1 12331 

31 

045 

855 

2METERS  ROAD 

31 

311.0 

SI 

00001  180 

112331 

31 

<J45 

305 

3'4ErE«S  ROAD 

31 

311.0 

SI 

00001 140 

112331 

31 

44 

305 

6C0UZENS  FrtY 

21 

21 

31 

311.1 

S 4 

00001200 

112331 

31 

124 

44 

IPURITAn 

21 

21 

31 

311.4 

SI 

0000  121  0 

1 12331 

31 

134 

124 

5L I V ERMJ I S 

21 

21 

31 

310.5 

SI 

00001220 

112331 

31 

1 1 1 

44 

50  .LODGE  FaY 

31 

31 

31 

ill  .0 

S 4 

00001230 

112331 

31 

2a  1 

1 1 1 

4J .LJOGE  FaY 

31 

31 

31 

sn.t 

S 4 

0000 1240 

1 12331 

31 

242 

281 

3J.LQCGE  Fay 

31 

31 

31 

310.6 

S 4 

00001250 

1 12331 

31 

2B5 

242 

20. LODGE  FnY 

4 

4 

31 

310.  5 

S 4 

00001280 

112331 

31 

285 

245 

1 GLEaOalE 

31 

31 

0.5 

SI 

00001270 

112331 

31 

245 

215 

2«04TEHEY 

31 

31 

0.7 

SI 

00001260 

112331 

31 

215 

175 

1 M(Jf<  TERE  Y 

31 

31 

0.4 

SI 

00001240 

1 12331 

31 

175 

185 

1DEXTE4 

31 

31 

0.5 

SI 

00001  500 

1 12331 

31 

185 

85 

1UUE4A  VISTA 

31 

31 

0 . 8 

SI 

00001310 

112331 

31 

85 

15 

1FULLEHTGN 

31 

31 

1.0 

SI 

00001320 

l 12331 

31 

15 

865 

2FULLE4T0N 

31 

31 

0.5 

SI 

00001 350 

112331 

31 

21 

855 

1SOOOLCRAFT 

31 

31 

31 

510.5 

SI 

00001 J40 

1 12331 

31 

41 

21 

13DAv I SO a AvE 

31 

31 

31 

311.2 

SI 

00001350 

112331 

31 

141 

41 

120 A v I SQM  AvE 

31 

31 

31 

510.5 

SI 

00001 560 

112331 

31 

221 

141 

llOAvISOf.  AvE 

31 

31 

31 

510.4 

SI 

00001570 

112331 

31 

251 

221 

100AV1SOM  AVE 

31 

31 

31 

310.5 

SI 

0000 1 580 

1 12331 

31 

242 

251 

40AVISON  AvE 

31 

31 

31 

310.3 

SI 

00001 340 

112331 

31 

31 

845 

5FE(.aELL  AvE 

31 

31 

31 

310.5 

SI 

00001U0O 

112331 

31 

101 

31 

OFEakELL  AvE 

31 

31 

31 

511.0 

SI 

00001«10 

112331 

31 

20  1 

101 

3FE.LAELL  AVE 

31 

31 

31 

510.3 

SI 

00001420 

1 12331 

31 

231 

201 

2FE.NAELL  AvE 

31 

31 

31 

310.3 

SI 

00001450 

112331 

31 

2el 

231 

IFEnkEll  avE 

31 

31 

31 

310.5 

St 

000014U0 

112331 

31 

21 

15 

l A YO’4  I iali  A VE 

31 

31 

31 

310.5 

SI 

0000  1 450- 

1 12331 

31 

31 

21 

2rt  YUM l vG  AvE 

31 

31 

31 

311  .0 

SI 

00001 4eO 

1 12331 

31 

44 

31 

SaYUMIKO  ave 

31 

31 

31 

310.1 

SI 

00001470 

1 12331 

31 

41 

85 

lLlVERMni3 

31 

31 

31 

310.  u 

SI 

00001480 

112331 

31 

101 

41 

2L I V ER  NO  I S 

31 

31 

31 

311.1 

SI 

0000 1440 

112331 

31 

1 1 1 

101 

3LIVE4N01S 

31 

31 

31 

310.1 

SI 

00001500 

112331 

31 

124 

1 1 l 

4L1v£rNUIS 

31 

31 

31 

310.4 

SI 

00001510 
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I 

I 


J 


1 12331 

31 

180 

191 

2DEXTEM 

31 

31 

31 

310.2 

SI 

00001520 

1 15531 

31 

191 

201 

30EX  TErf 

31 

31 

31 

310.9 

SI 

00001530 

112331 

31 

215 

221 

IlI*«OOD 

31 

31 

31 

310.6 

Si 

00001540 

1 12331 

31 

221 

231 

2L  IN-'iOOU 

31 

31 

31 

310.7 

SI 

00001550 

112331 

31 

245 

251 

1 12TH 

31 

31 

31 

310.6 

SI 

0000l5o0 

112331 

31 

251 

261 

2127* 

31 

31 

31 

310,5 

SI 

00001570 

999 

00001580 

111176 

76 

383 

153 

7E IGHT  «I 

RD 

7a 

76 

3 

30.8 

SI 

00001600 

111176 

76 

783 

363 

b£ I Gh  T * I 

RO 

7 6 

76 

3 

30.8 

Si 

00001610 

111176 

76 

783 

765 

3FI4ELANE 

89 

89 

76 

762.2 

Fl 

00001620 

111176 

76 

765 

775 

47EN  mjl£ 

90 

76 

760.4 

Si 

00001630 

111176 

76 

775 

165 

IFlNELAl.e 

76 

7b0 , 9 

Fl 

00001640 

111176 

76 

165 

153 

8L I VESNO I S 

76 

76 

76 

76  1.6 

SI 

0000 1 o50 

111176 

76 

383 

le5 

5yoOD.*arO 

Av 

76 

76 

76 

761.7 

St 

00001 660 

999 

0000lb70 

111169 

89 

au3 

783 

5EIGH7  Ml 

80 

89 

89 

3 

30.8 

Si 

OOOOleSO 

111189 

89 

555 

443 

4tIGM7  i*  I 

90 

69 

89 

3 

30.a 

SI 

00001690 

1 1 1 1S9 

89 

565 

555 

SOEmU  I filORE 

8 

89 

89 

2.0 

SI 

00001 700 

111189 

89 

5o5 

475 

1 t E 4 mile 

90 

89 

891.0 

SI 

00001710 

111169 

89 

475 

465 

2TEi4  mile 

«0 

89 

890.4 

SI 

00001720 

111189 

89 

465 

7oS 

3TE  J mile 

90 

89 

890.1 

SI 

00001730 

111169 

89 

443 

451 

1C*MYSL£9 

F Y 

39 

89 

'89 

391.1 

$4 

00001 740 

111189 

89 

451 

«e8 

2S  r£R*E^«SGN 

89 

89 

84 

841.3 

3 1 

00001750 

111189 

89 

475 

45  1 

1J0M.4  8.  90 

89 

89 

89 

891.2 

SI 

00001 760 

499  0000  1 7 70 
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I 

i 

i 

l 

iO  0 0 40  0 0 140  0 0 210  0 0 240  0 0 ilO  0 0 760  0 0 0000001O 

840  0 0 *1  00000020 


4-3  Listing  of  Records  in  Input  Data  File  IN,  Control  Card  File 

Identified  In  JCL  as  FT01F001.  Detroit  Test  Case,  using  all  unit 
areas  In  Links  File 
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Section  V:  Operating  Instructions 

Since  the  computer  installation  at  which  the  Transportation  Submodel 
is  to  be  used  is  different  from  the  one  at  which  it  was  developed,  only 
general  guidelines  are  offered  on  how  to  assemble  and  use  the  FORTRAN  source 
cards.  Specific  local  operating  instructions  can  be  developed  by  the  local 
programming  support  staff. 

The  FORTRAN  source  cards  for  the  Transportation  Submodel  comprise  5 
distinct  sets  of  main  programs  with  associated  subroutines.  The  program 
names  are  NETWORK,  ALLNET,  DAMAGE,  FLOYD,  AND  TVLREC.  Three  different 
methods  of  assembling  and  using  the  complete  package  have  been  used.  One 
method  is  to  create  a separate  load  module  for  each  of  the  5 sets  and  to 
execute  them  in  succession  as  a series  of  steps  in  a job.  The  job  control 
language  (JCL)  for  executing  the  load  modules  is  shown  in  Figure  5-1.  In 
creating  the  load  modules  provision  must  be  made  to  provide  the  multiple- 
use  subroutines,  NETPRN,  STORE,  and  UNMASK  either  as  load  modules  in  a job 
library  or  by  creating  duplicate  source  decks. 

In  the  second  manner  of  assembling  the  programs,  ALLNET,  DAMAGE,  FLOYD, 
and  TVLREC  are  compiled  as  subroutines  called  by  NETWORK.  In  this  case  the 
necessary  statements  CALL  ALLNET,  CALL  DAMAGE,  etc  are  added  after  the  state- 
ment 30  CONTINUE  in  NETWORK.  The  statements  SUBROUTINE  ALLNET,  SUBROUTINE 
DAMAGE,  etc.  are  added  to  the  source  decks  and  STOP'S  are  replaced  by  RETURN' 
No  parameter  list  is  necessary  in  the  calling  or  subroutine  statements.  The 
JCL  for  a compi 1 e-1 i nk -edi t-go  in  this  situation  is  shown  in  Figure  5-2. 

The  third  manner  of  assembling  the  programs  is  a variation  on  the  second 
described  above.  NETWORK  is  converted  to  a subroutine  in  the  same  manner  as 
the  other  programs.  The  whole  package  of  subroutines  is  called  by  a simply- 
written  main  program.  The  primary  feature  of  the  main  program  is  that  it 
contains  a loop  to  simulate  time  iterations  where  the  Resource  File  is  updated 
with  different  values  of  BOS.  This  third  technique  of  program  assembly  is 
mentioned  without  an  example  just  to  suggest  what  might  be  done.  In  fact,  it 
was  done  with  the  programs  in  an  earlier  stage  of  program  development  where 
the  Master  Status  File  was  used  instead  of  the  Resource  File. 


//PROCESS  JOB  RTI ,C43.P3495e, J«D,m=1 , PWT Y=1 , T = 1 , F QRmSs’ah I TE 

//•  RTI.C43.PC495o.JrtO.  PROCESS.  CNTl  TO  EXECUTE  l 9 A NSPOfi  T A T I UN  SUBMODEL 

//•  4 S A SEQUENCE  OF  LOAD  MODULES  IM  JOB  STEPS. 

//SI  ExEC  PGM='NETrtORK,REGION=9BK 

//STEPLIB  OD  DISPsSriR,OSNs»U.CU3.P049S6.J*D.NETrtgR*.LL)AD 
//FT01F001  DO  0I3P=SMR,0SN=RTI .C43.P049S6. JAO.F TO1FO01.CATA, 

//  DC9  = BuFvi)=l 

//FTOJFOUl  QD  SYSUJT=A,DC3=a  JFnO=1 

//FT30F001  DO  OSNsRT I.C43.PG49Bo.JrtD.NErH3T.0ATA,0 ISPs (NErt.C A TLG)  , 

//  UNIT  = 01  Six.  VUL=3ER=RT  1222,  SPACE *(  TRk,  C 1 , 1 } ,RLSE)  , 

//  DCB=C9tCFR=VBS,LSECL=lb08,6LiX3lZ£  = 1612,dUF-',O=n 

//FT53FOOI  DU  DISP  = (NE.i.C4TLG) ,CSN=RT I .Ca3.PU495o.J, vD. PATHS. DATA, 

//  UNlTsDISiX,  VOL  = S£H  =R  T 1 222,  SPACE  = (TRK,  (50,50) , RUSE)  , 

//  DCBs  CRECFMsvBS.LRECLsl  3022,  BLXSIZE  = 1 3026, BUFivOs  i ) 

//FT51F001  DD  0 ISP=SHW , OSNsR  r I .C 4 3 . P04956 . >0 ,L INKS . OA T A , DC B=BUF xC  = 1 
//S2  EXEC  PG*  = ALL,iET,REGIQN  = 2S2K 

//STEPLIB  DD  OISPsSHR.DSNsHTI  .Cu3.PDU95o.  JiiD.ALL  ,ET  .LDaO 
//FT03F001  DO  S r SO JT s A , DCB=8 JFNOs 1 

//FT30F00  1 DO  DISPsSHR,DSN=Hri.C4  3.P0U95o,JviD..\£TLIST.DArA,0CB  = euF  xU  = 1 
//F  T5  3F00  1 OD  0 ISP  = SHfi , OSNsR  T I . C 4 3 . P'J495o  . J.xO  . P A T hS  .DA T A , uCBsfluF mO=  1 
//FT70F0UI  DO  0SN=ii42B0UN, UNIT=0ISK,DI3P=(, DELETE) , 

//  SPACEs( TSr, (50, 1 ) ,RLSE) , 

//  OCBs (RECF 4 = VBS,LRECL=lb08,BLtxSIZE=l  1260) 

//FT71F001  DO  0SNsiitt2B0UN,UMTs0ISix,0ISP=(, DELETE), 

//  SPACES  (TRix,  (50, 1)  ,RLSE)  , OCBs  ».FT  TOF  (10  1 
//FT80F001  OD  0SNS4480UN2 A , UN  I TsD I SK , D I 3?  = ( , DELE TE  ) , 

//  SPACES  ( TRIX , (50, 1 ) , RLSE)  ,DCB  = ».F  T70F001 
//FT81F001  DD  DSNs A&3UUN26 , UN  I T = 0 1 SK , 0 1 SP= ( , DELE TE ) , 

//  SPACE=(TRix,  (50,  l)  ,RLSE)  ,DCB  = » .F  T 70F00  1 

//FT91F00  1 DD  0 I SP  = ( VE* , C A TLG ) , OSNsRT  I . C«  3 . P04956  . Ji.D  . A280UN  .DA  T A , 

//  SPACES (TRK, (SO,  1 ) , RLSE ) , uN I TsO I SK , VOL  = SE R = R T 1 222 , DC Bs «. F T 7 OFO 0 l 
//F  T92P0OI  00  DISP=(Nt«,C4TLG)  , US N=»  T I . C 4 3 . P0«95o  . Ji.O  . 80  Jn28  . 0 A T 4 , 

//  SPACES  ( TRK  , (50, 1 ) ,RLSE)  ,OMT  = DISK,  VUL  = SER  = R T 1 222 , 0Cb  = • . F T 70F00 1 
//S3  ExEC  PG-isOAMAGE , hEGIQN=80iX 

//STEPLIB  DD  DISP=SHR,DS'4sPT I .C43.P0U956. J«D.DAh4GE.L3AD 
//FT03F001  OD  SYSQJTsa.DCosbJFNO*! 

//FT30F001  DO  OISPsSHP, OSNsRT  I . C 4 3 . P0495o  . J*0 . f.E  T L I S T . D A T A , OCB  sriJF  N0  = 1 
//FT91F001  DD  DS  <*RT l .C4J .P04956. Jfl0.A260UN.DAT A , DISFsSH* , DCB  = «UF NOs  l 
/ /r  T92F  00  1 DD  DSnsRI  I . C4 3 . PU495o . J *D . B0UN2B . 0 A T A , 0 1 SPsSHR , OCBsBUF  ,y=  t 
//A  T93F001  OD  0I3PS(  NE/i.CATLG)  .OSNsRT  I . C 4 3 . PQ4956  . J.'iC  . D I 5 T ANCE  . 0 A T A , 

//  SPACES ( TRK, (50, 1 ) , RLSE) ,UNl TsQISK, VUL=SER=RT I222,CCB=* .FTPlr 001 
//FTblFOOl  DO  0 I SPsjHR,  DSnsR  T I .L  43. PU49S6.J.-.0.  rESOUhCE.DAT  A2, 

//  OCBsBUFNJsl 

//S4  EXEC  PG'-'sFL J T D,REGlON=100ix 

//STEPLIB  OD  OISPsSHR.DSNsRTI .C45.P0«9So.JrtD. FLOYD. LOAO 
//FT03F001  DD  3 Y 30 J T= A , UCB=BUF NO = 1 

//FT 30F00 1 JO  oISPsSHR, OSNsRT I . C 4 3 . P04956 . J rtO . NE T L I S T . D A T A , DCs seuF N0=1 
//F  T9  3F00 1 DO  OSNSRT I . C4 3 , P04956 . JnO , DISTAnCE.DAT A, 01 SP=SriR , DC d sBuFNU* 1 
//FT55F001  DO  DISPs(N£rt,CATLG), DSN sWTI,C43.P049S6,J'rtD, FINAL. DATA, 

//  SPACES (TRK, (SO, 1 ) ,RLSE) ,UM  TsDISK, vOL*SER  = R T I 222 , DCB=» .F T93F 00 1 

5-1  Llstlnq  of  JCl  to  Execute  Transportation  Submodel 
as  a Sequence  of  Load  Modules  In  Job  Steps 


oooooo i o 
00000020 
00000030 
00000040 
00000050 
00000060 
00000070 

ooooooao 

00000090 
00000100 
00000110 
00000120 
000001 30 
00000140 
00000150 
00000160 
00000170 
00000180 
00000190 
00000200 
00000210 
00000220 
00000230 
00000240 
00000250 

0000026a 

00000270 
00000290 
00000290 
00000300 
00000310 
00000320 
00000330 
00000340 
00000350 
00000360 
00000370 
00000580 
00000390 
00000400 
00000410 
0 U 0 0 0 4 2 0 
00000430 
00000440 
00000450 
0 0000430 
00000470 
00000440 
00000490 
00000500 
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//FTbOFOOl  00  OSN=&i'*£XUS£,UNITsJIS*,OISP*(,OfcLfcrE), 

//  SPACEs(TWK, ( 1 , 1 ) ,HlSE)  , 

//  OC3=(PECFMsvBS,LRECLsieOfl,3l.4SlZE=le!2) 

//FT7OF00  1 00  OS.'.s&iTE  XP  t , dN  l T SO  I S«v  , J I SP=  l , dELE  T E ) , 

//  SPACE=irMK,(50,l),RLSE), 

//  OCB=f «ECFMsvSS,LHtCLslo08.aLK5IZt=l  1E60) 

//FT30F001  00  0SN=44 T E*P2 , UN  I T so  I S* , 0 1 SPs ( , OEL t T E ) , 

//  SPACE  = crp*((5o,u,PLSE),oces«.Fr7oFooi 
//35  ExEC  PO.-'srvLHEC.PE&IONsloOK 

//STEPLIB  00  JlSPsSMR, OSMsK  r I . C«S . P3<*956.  J*D . T VL7EC  . LOAO 
//FT03F0Q1  00  S»3OjT=4,0Ca=ajF'*0sl 
//FT30F001  00  0S'-s4&HETI0.OISPs(3L0,OELETE) 

/ /F  r 30F 00  1 00  OISPsSHR.OS'iaHII  .C  J3.P04956.J.-.D.. .ETUIS  7.04  r a,  CCBsBUFnCs 
//FT55F001  00  0ISPsSMH,0SN=RTI.Ct*3.P0«956.jAO.H'l4L.04T4,0CB=aoFNdal 
//F  T9«F  001  00  0 ISP*  C '*£<’> #C4TLG)  .OS^aRT  I ,C«3  ,PU«95b . J»0 . T vLREC  .OAT  4 , 

//  SP4C£  = (TH* ,(50,1 ) ,«LSE)  ,UN1  TsOlSix,  VUL  = SERs«ri222, 

//  0C8=(RECFMsFe,LRECL=20,aL4S I ZE® 130201 
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00000510 
00000520 
00000530 
OOOOOSjO 
00000550 
000005n0 
00000570 
00000530 
00000570 
OOOOObOO 
OOOOOol 0 
00000620 
1 00000630 

000006 JO 
OOOOOeSO 
0000066 0 
0000Ob70 


j. 


I 

I 


I 

I 


//EXECUTE  JUd  pTl.C4i.P0495d,J/.0,-<=l,PRTY  = t,r:2 

//*  RU. C4i.PQ4956.JwQ. EXECUTE. C *Tl  TO  EXECUTE  T 9 AN SPQP T * T 1 UN  SUBMODEL 
//*  AS  A SEQUENCE  UF  calls  «ITnlN  'NET.-.PK'  To  Th£  OThE*  PROGRAMS  AS 
//»  SUBROUTINES 
//SI  ExEC  PThC 

//C.SySIN  00  OSt.sHTI  ,CUi.P0495o.  Jw0..iETn0«K  .POP  T , 0 1 3P=Shk 
//S2  ExtC  p T hC 

//C.SYSIN  00  QSN  = PTI.C43.?U495b.Jp0.ALLNET.FURT,:)I3P  = 5,-w 
//S3  txEC  PTHC 

//C.SYSI*  00  OSNsRTI.C43.PO«956.J«O.OAMAG£.rURT,DISP=SHP 
//S4  EXEC  PThC 

//C .SYS IN  DO  DS  cPTI  .C4i.P0495b.J, *0. FLOYD. Fl)HT,DISP=SH9 
//S5  tx£C  FThCLG,R.G=500p 

//C.SYSIn  00  0Sfts»TI.C4i.PU495o.Jrt0.TVLREC.FQRT»5ISP=SnP 
//G.FT01F001  00  OISP  = SH«,OS »=S TI ,C4i.P0495o. JmO.PTO 1P00 1 .DATA, 

//  OCB*QjF*U=l 

//FTOiFOOl  00  SYS0UT=A,DCBs8UFNC*l 

//FTJOFOOl  00  0SNsRTI.Cu3.PO995b.J*D.N£TLIST.DA7A,DISP=C  .Ey.CATLG), 

//  UMTsDISK,  VUL=SER=P TI 222, SPACE =(TRK  , ( 1 , 1)  ,AlSE)  , 

//  0C5=(nECF.4=VBS,LRECL=lbC8,BLPSIZE  = lot2,dUFN0=l  ) 

//FT53FQ01  00  0 ISP* ( Ntw, C A TLG ) , QSNs« T I . C 4 3 . PU495b . Jv.D  . P A THS ,0A T A , 

//  UNI 7=0 IoK,VUL=SEP*R 71 222 »SPACE  = (TkK, (50,50 ),«LS£) , 

//  DC»=(»ECF.4=vtJS,LRECL=13022,dL<SIZE  = 130  2b,auP''.0=l  ) 

//FT51F00  1 00  JISPsSHH,DSN=«TI  ,C 4 3 . PJ4956  , JwO  ,L I VPS  . OA  T A , OCB=diJF N0  = l 
//FT70F001  00  OS *=44A2BUUN,UNl TsOISK,OlSP=( .uElETE) , 

//  SPACE  = URK,  (50. 1)  ,PLSE) , 

//  OCB=CRtCF  ' = VdS»LPECL=l608,dLPSIZE=l  1260) 

//FT71F001  00  0SNsi&b2di)UN,UMITs0ISix,0ISPs(, DELETE), 

//  SPACE: C TP*. (50, 1 ) ,HLSE) ,OCB=».FT7QF001 
//FT30F001  Ow  OSN  = 4AdUUN2A,UNIT  = OISK,OISP=(,OELE'TE), 

//  SPACEs(TPP,(50,t),PLSE),OCa:»,FT70F001 
//FT81F001  00  0SN  = 44d0UN2B, UNI T = JISK ,01SP=( .DELETE)  , 

//  SPACE® (TAP, (50, 1) ,HLSE) , 0C3: » .F T70F00 t 

//FT91F001  00  OISP:(',£a,CATLG) ,OSn:PT I . C « 3 . P04956 . J*C . A2dUUN.DAT A, 

//  SPACE*(TWK,(50,l) , PLSt ) .ON  I 7 = 9 1 SK , VUL=SEr=B T 1 222 , OC 5 = * . F T 70P00 1 
//FT92F00  1 00  0 1 SP  = ( NEW , C A TLG ) , OS*=B T I . C4 3 . P(J495b . JnO . dOjw2a . 0 A I A , 

//  SPACE=(  TP*,  (50, 1 ) ,9LSE)  ,.UW  I T:0 1 SK , vQL=SER=h  T 1 222 , 0C3  = * . F T 70F  00 1 
//FT9  3F001  00  JISP  = ( nE.v.CATLG) ,OSN  = HTI .C43.PQ4956. JpO.OISTAr  Ct.UAfA, 

//  SPACE: ( TPK, (50, U ,«LSE ) ,UN I 7 = 01 SP , VOL *SE * =9  7 1 222 , DC c = * . F 7 9 1 K 0 0 1 
//F761F001  00  01SP:SnP,OSi'J:PTI.C4i.PQ4956.J’Y0.PtSJU9LE.OA7A2, 

//  OC b:0UFNL's  1 

//FT55FO01  00  .0ISPs(NEw,CA7LG) ,0SN*P 71. C43.P04956.  J«0.FI  iAl.OATA, 

//  SPACE: ( 79P, (50, l) ,HLSE) ,UNI 7:0 1 SP , VOL»SE«s« T 1 222 , 0Cfl  = * . F T93F00 1 
//FTbOFOOl  00  DSN  = 44  w EXUSE, Ui.IT  = 9ISP,DIS?  = (»UtL£TE), 

//  SPACE:( 7HK ,(1,1 ) , PLSE ) , 

//  0C8=(P£CFM:vdS,LPECL=lb08,6LPSIZE:lbl2) 

//  SPACE: (TPK, (50, 1 ) , PLSE ) ,0C8: * . F T70F00 1 

//FT94F001  00  0 1 SP«(Nfcw,C A TLG ) ,OSN*RT I.CaJ.PU4956.JwO.TvLMEC.DAT A, 

//  SPACE* (TPK, (50, 1) .RUSE) ,U  UTsOISp. VQL:SEp:PTI222. 

//  DCd=  ( PECF'**F  d,  LRECL  *20, 3LP  SIZE*  13020) 


000000  1 0 
00000020 
00000030 
00000040 
00000050 
00000060 
00000070 
00000030 
00000090 
OOOOOIOO 
000001  10 
00000120 
00000130 
00000140 
00000150 
00000  tbO 
000001 70 
000001 80 
00000190 
00000200 
00000210 
00000220 
00000230 
000002*0 
00000250 
00000260 
00000270 
0000-0230 
00000290 
00000300 
00000310 
00000320 
OOOOv  330 
000003*0 
00000350 
000u03b0 
00000570 
10000330 
00000390 
00Q0O400 
00000410 
00000420 
00000430 
00000440 
00000450 
000004b 0 
00000470 
00000*30 
00000490 
00000500 
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Section  VI:  Suggestions,  Warnings,  and  Changes 

Several  improvements  could  be  added.  At  one  stage  of  development  the 
user  had  to  supply  the  list  of  nodes  connecting  unit  areas.  In  an  effort 
to  reduce  the  number  of  files  and  the  amount  of  manual  data  preparation 
required,  it  was  decided  to  let  the  computer  do  the  work  of  determining  unit 
area  connectivity,  and  this  file  was  eliminated.  As  experience  suggests  and 
familiarity  with  the  program  and  files  increases,  it  may  become  desirable  to 
reinstate  this  list  to  provide  an  override  option  to  describe  linkage  between 
unit.areas. 

Another  improvement  would  be  to  permit  changing,  say,  road  type  from 
major  freeway  both  ways  to  one-way  for  "movement  to  shelter"  procedures  during 
the  course  of  the  time  iterative  scenario.  At  present,  the  Links  File  which 
describes  the  network  in  its  initial  state  is  "outside"  the  time  iterative 
loop.  Thus,  the  only  way  to  effect  a change  in  states,  other  than  damage,  is 
to  so  describe  the  network  initially  in  the  Links  File. 

Another  improvement  would  be  to  create  a second  version  of  the  step 
that  computes  the  shortest  paths  between  unit  areas  when  the  total  number  of 
unit  areas  is  less  than  approximately  150.  This  second  version  would  replace 
FLOYD,  which  uses  the  distance  matrix  stored  on  disk  or  tape,  wi;h  something 
similar  to  SPATH,  which  stores  the  entire  matrix  internally  in  core,  at  a 
great  savings  in  I/O  time.  Hand  calculations  suggest  that  150  unit  areas 
will  generate  distance  matrices  (150  x 150)  just  within  core  limitations. 

This  improvement  would  be  justified  if  there  were  anticipated  a large  number 
of  study  areas  containing  less  than  150  unit  areas.  Another  advantage  to  be 
gained  would  be  that  Dijkstra's  algorithm  could  then  be  added  to  resolve 
for  oriqin/destination-specific  shortest  paths  after  damage. 
t Some  infrequents  used  programming  techniques  are: 

. (1)  the  use  of  9999  to  represent  a distance  of  "infinity"  between  nodes 

(2)  storing  multiple  values  in  pre-defined  bit-fields  of  a single  word 
(see  the  description  of  STORE  and  UNMASK). 

(3)  storing  the  value  of  distance  in  the  upper  half  of  a word  and 
policy  in  the  lower  half  for  words  in  the  distance  matrix.  The 
value  216  is  used  arithmetically  to  shift  and  extract  values,  and 
will  be  seen  to  occur  frequently  in  the  program. 


A technique  for  changing  recursively  between  two  device  numbers  in 
READ  and  WRITE  statements  is  used  in  FLOYD  and  ALLNET.  In  order  to  orocess 
data  written  to  an  output  device,  (e.g.,  number  J),  the  device  number  in  the 
READ  statement  must  be  set  to  value  J,  and  the  new  WRITE  output  device 

number  is  set  to  value  I,  the  prior  value  of  the  input  device.  Thus,  if  we 

are  reading  from  I and  writing  to  J,  the  recursion  is: 

(I  + J)  - I ^ I 

(I  + J)  - J > J 

i 

The  most  likely  error  to  occur  will  be  an  inacessible  node  in  a unit 

area  network.  This  condition  will  be  noted  by  the  value  9999  appearing  some- 

where in  the  shortest  paths  matrix.  Examine  the  Links  File  for  that  unit 
area  for  transposed  digits  in  link  and/or  node  identification  numbers,  an 
incorrect  node  i.d.,  inclusion  of  a node  from  another  network,  or  incorrect 
specification  of  a one-way  street. 


Program 

Listing 


i 


RTI.C4  3.P0u<J5b.JwD.NET..0HK.F0RT 

PROGRAM  NETWORK  In  TRANSPORTATION  SUBmQUEL. 


«E^OV£  *c*  IN  SUBROUTINE  NETwRK  TO  USE  AS  A SUBROUTINE.  ALSO  CHA  ,GE 
sro?  ANU  RETURN  CASUS  BELOw. 

SUBROUTINE  NETwRK 
CQMMUn/IO/IN, IUUT 

COm“UN/olOCM/NINE  le,  INF  I N , K I NF  I N , ITwUlb.NETSlROO) 

„ ' COMMON/6UOCK2/OI3TAN(  75, 75 ) , IK/LFL ( 75 , 75 ) , NO  THAN ( 7b J , I\DE *T ( 75) 

integer  hiahchis),  oistan 
I T *0 1 o a 2**10 
NlNElb  a BBRB 
INF  IN  a NINE lo 

KlNFIN  a inf  n»  IT  *010 

HEAD  IN  LIST  OF  NETWORKS  TO  BE  PROCESSED. 

CALL  CNTLHD CNUMNET) 

AT  TbIS  POINT,  wE  HAVE  ALL  NtCESSAHT  DATA  To  BEGIN  Th£  SHORTEST  path 
PROCEDURES.  ThET  »IlL  o£  IMBEDDED  wIThIN  a LOOP,  SO  That  EACh 

network  will  be  analyzed  in  turn. 

DO  10  IA1,NUMNET 

CREATE  MATRIX  OF  DISTANCES  BETwEEN  DIRECTLT  LINKED  NODES  IN  The  I-Th 
NE  T wORk  • 

CALL  SETuPCI, NODES, 450) 

PRINT  The  INITIAL  POLICY  MATRIX  BETWEEN  NODES  FDW  I-TH  „ET.-.()RK  ON 
DE  VICE  ' IOUT ' . 

CALL  PRNOUTCIPOLFL, NODES, I ) 

PRINT  The  INITIAL  MATRIX  OF  DISTANCES  BET..EEH  NODES  FOR  THE 
I-Th  NETwORk  On  OtViLE  ’IOuT' 

CALL  PRnOuT (OISTAN, NUDES, I ) 

COMPUTE  The  paths  of  MINIMUM  DISTANCE  BE  T„EEn  all  ,uLES  IN  The  I-Th 
NE  T WORK  : 

CALL  SPATHC  NODES) 

SAVE  the  RESULTS  OF  The  SHORTEST  PAThS  COMPUTATION  UN  FILE  5J  FUR 
SUBSEQUENT  PROCESSING  OF  SHORTEST  PATHS  BE  T wt£N  UnIT  AREAS. 

7-1  Listing  of  Fortran  Source  Program  NETWORK 
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/•airE(5i)Ners<i)*nooes,*orwA>M,DisTAn 

c 

C PRINT  T«fc  FINAL  SOLUTION  POLICY  A„D  DISTANCt  MATRICES  3N  CEvICE 

c 'lour', 
c 

call  prnouuipqlfl. nodes.  n 

CALL  PRNOUT(OISTAN,NOOES,I) 

30  CUNTInoE 
REaINO  S3 

c 

C have  PROCESSES  NJMNET  Jfiir  area  NETmORhS . SAVE  THE  P A C a E 0 nETiGRa 
C I.O.  NUMBERS,  IN  NErS(I),Ist,2,...,NU«NET,  UN  FILE  3C  . 

c 

flR I TE ( 30 ) 4UMNE  T ,N£  TS 
RE "I MO  30 
C 

C REMOVE  'C*  IN  COLUMN  1 TQ  ExECuTF  REST  OF  T ft  ANSP  OR  T A T I UN  MODEL 
C PROGRAMS  AS  SUBROUTINES. 

C 

C CALL  AlLNET 

C CALL  OAMAGE 

C CALL  FLOTO 

C CALL  TvLREC 

STOP 

TO  EXECUTE  '..ETiORK*  AS  A SUBROUTINE  remove  'C*  In  RETURN  CARO, 
REMOVE  'STOP*.  AND  REMOVE  'C'  In  SUBROUTINE  .ET.-.QRa  CARO  ABUvE. 

RETURN 
EiVO 

block  oata 

COMMON/IO/IN, IOUT 
DATA  IN/I/, IJUT/3/ 

END 

subroutine  se tup (I  am, nodes, «) 

C........ ................................. ............................. 

c this  subroutine  reads  the  network  data  for  netmjm  am,  creates  the 
c INTIAL  matrix  OF  distances  OETflEN  DIRECTLY  lINaED  NODES. 

C IT  ALSO  PACKS  The  ZONE,  EOC,  GROUP,  SECTOR,  and  unit  area 

C I.O.  NUMBERS  INTO  .ETSCIAM),  *hiCh  EnTERS  SETJP  containing  Only 

C the  NEThORA,  £.G.  unit  area,  I.O,  NUMBER. 

C VARIABLES: 

C NO TRAN 1 1 ) : ACTUAL  NOOE  NUMBER  (ThE  SHORTEST  path  ROUTINES  •'ILL  aGRa 
C «ITH  dummy  NUMBERS  l,  2,  ETC.) 

C OISTA.N(I,J):  distance  from  node  I TO  node  j. 

C IPQLFL  :„OR*ING  STORAGE  FOR  SURT;  BECUMES  IP0LFL(1,J)  In  'SPATh', 
C : Th£  POL  ICY  MATRIX. 

C INOExT  5 aORkING  STORAGE  FOR  SORT. 




0U000510 
00000524 
00000530 
00O00340 
00000553 
300 005a 3 
00000570 
00000580 
000Q05RO 
OOOOOaOQ 

oooooa  i o 
00000620 
00000630 
00000640 

oooooaSo 
00000660 
03000670 
OOOOOoBO 
00000690 
00000  700 
00000  71  o 
00000720 
00000730 
00000740 
00000750 
00000760 
00000770 
00000730 
00000700 
00000800 
00000310 
00000820 
00000830 
JOOOOauO 
00000350 
0 0 y u 0 3 3 J 
00000370 
30000830 
00G0O8R0 
0000090  3 
0000091 3 
00000920 
00000930 
0oQ0Q9u0 
00000930 
0 0 0 0 0 9 3 0 
0uo0097o 
00000930 
000009RO 
00001000 
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COMMON/ IQ/IN, lOUT  00001010 

COMMON/bLUCR  t/NtNEl6,INFIN,KINFI,N,  I T rtU 1 6 . NE T S ( 40 0 ) 00001020 

C0MM0N/dLUCK2/0ISrANC75, 75) . IPOLFU7S,  75)  . INQE  * T ( 75)  . LOTS  AN(  75  ) 00  00  10  JO 

DIMENSION  InDU£S(75),hIARCh(5)  00001040 

INTEGER  FOU  NO.  »E*  I NO#  ZONE  I D . EOC 1 D , GRUP I D , SEC T 10 , H I A RCM, 0 1ST  AN  00001050 

0000 1 060 
00001070 

IONE  - 1 00001060 

.ETNUM  = NETS(IAM)  0009 1 ORO 

00001 100 

INITIALIZE  THE  DISTANCE  MATRIX.  00001110 

OOOOl 120 

DO  15  1=1,75  00001130 

DU  14  Jsl.75  00001140 

IPOLFL(I.J)  = NINEla  00001150 

14  DlSTAN(I,J)=NlNEle  00001160 

IPOLFLd.I)  * 0 0000  1 1 70 

15  DISTANd.i)  s 0 0000  1 130 

00001 1 RO 

INITIALIZE  the  matrix  of  NODE  I.J.'S  00001200 

DO  lo  LSI, 75  00001210 

I NGOES CL  1 = L 00001220 

INOEXT(l)  = 0 0000  1250 

la  NOTRAN(L)  s 0 00001240 

NODES* l 00001250 

e5  FQUNOSO  00001260 

REaK.O  = 0 00001270 

00001230 

read  the  LINKS  file  - LOOKING  FOR  nETaCRk  I.D.  NUMBER  nETSCIAM)  00001240 

00001 300 

70  REAO(51,80,ENOs71  )hIARCH,IFOH»n,L,LF  DR , lb AC* , L VLBAC , L 1 NkNU , Am ILES . 0000  1 3 1 0 

1 T y PE , I T rPt  0000  1 320 

30  F (JR MATC2I1, 212/13,2(13, II), l4,2JxFi. 1,1* A 1/11)  00001  330 

NET  NO  * HIAHCHC5)  00001  540 

IF (NETNO.E3.NE TNUM)GU  TO  31  00001354 

IF ( n£ TNC.E3 .RRR . AND.FUUnO . E 3 . 1 ) GU  TO  Uo  00001560 

GO  TO  70  00001 570 


71  RE.nINO  51 

IF (FOUND. £3.1) GO  To  lie 
RE*  I NO  * RE..INO*l 
IF  (RE.nINO. LT.2)G0  TO  70 
rR1TE(I OUT, 33) nETNom.nET NOR 
RETURN  l 


OOOOl  320 
00001 330 
00001 340 

0000  1 354 
0000 1 560 
00001 570 
00001 580 

00001  JRO 
00001400 

0000  1 4 1 0 
4000 1420 
00001430 


38  FORMAT!/*  IN  SobR.  SEToP.  2.0  REwInO  ON  LINK  FILE  (DExiCE  51)  L00u000l440 


IKING  FOR  NETWORK  NUMBER ',14,' 


NOT  FOUND 


1HOOT  COMPUTING  SHORTEST  PA ThS  FOR  NE  T..U«k  * , I 4 ) 


81  FOUND* l 

CALL  MA  TCH(  IFOR*,  lOfxE,  NO  TRAN,  NOOES,  IX,  I ,a,ATCh  1 ) 


RETURNED  TO  -'AI,  *1700001450 
I 00001460 

00001470 
OOOOl  490 
OOOOIRPO 

1)  00001500 
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nonnnnnon 


c 

c 


c 

c 


iF(M4rcMj.£u.nGa  to  oa 

NDTRANIhOOES)  = IFORW 

i * nodes 

NOOESsnQCESn  l 

aa  call  match(I6ack» iaf.e,NorsA.M,. mooes, I*. J^-’ATcwa) 


IF ( MA  rcn2 . ED • 1 )GQ  TO  S3 
VOTRANCNGDES)  = ISACk 
JsMOOtS 
nuDES=nODES»1 

68  IMILES  = A«ILES*10.  t 0.5 

IF( (ITrPE.EQ.3) .UR. dTYPE.E3.6)  ) GO  TO  110 
90  OISTAN(I.J)  = IMILES 

IF(ITYPE.tU.2.0R.ITYPE.EQ.5)GO  73  70 
110  0ISTAN(J, I ) = IMILES 
GO  TO  70 


116  NODES  = NODES” 1 
OO  10  1=1,5 

10  NETSCIAX)  = SrOHE(HlA«ChCI),NErS(IAM),I j 

call  rank (ione, nooes, nOTRAn, Index  t , I nodes) 

03  40  1=1, NODES 
NGTRAN(I)  = INOExI(I) 

ILsINOOESCIJ 
oo  <*a  j=i,nooes 

JL  = INOOESCJ) 

«0  IPOLFLtl.J)  = OISTANIIL, JL) 

OO  41  1=1, NODES 
00  41  J=l,NOOES 
OISTANd.J)  = IPOLFUI.J) 

IPOLFUI.J)  = 0 

IF  (OISTANd.J)  .NE.NInE16)IP0LFHI,J)=J 
41  CONTINUE 
RE  TUHN 
END 

SUBROUTINE  SPA 7m ( yQOES ) 


THIS  SUBROUTINE  JSES  FLOYO'S  AlGJR I Thm  TU  F I i,o  TmE  SHORTEST  paths 
FROM  all  NODES  TO  ALL  OTHER  NUuES  (OH  FROM  ALL  NODES  TO  A 
SPECIFIC  DESTINATION  IF  JsCONSTAnT.) 

SEE  'A,.  APPRAISAL  OF  SOME  SHOW  TEST-PATH  ALGORITHMS'  oi  STUART  E 
OREYFUS,  0RSA-l9O9,  MAY-JJNE,  VOLUME  17,  NO-SER  3,  PC.  396-413 
VARIABLES  : 

UISTAN  S DISTANCE  MATRIX  FOR  FlOYD'S  ALGORITHM 


implicit  integer  t o > 

COMMON/ 10/ IN, I OUT 

COMMQN/BLQCK  l /NINE  16,  INF  IN,  K IMF  IN,  I T.<ulb,.,£  TSt  40  0) 


00001510 
00001520 
00001530 
00001540 
00001550 
00001560 
00001570 
00001580 
00001590 
00001600 
OOOOlolO 
00001620 
OOOOleSO 
0000 1640 
00001650 
0000 1 ooO 
00001670 
00001630 
00001 o90 
00001700 
00001710 
oaooi 720 
00001 730 
00001 740 
00001 760 
00001 760 
00001770 
00001780 
00001790 
00001800 
0 0 0 0 1 3 1 0 

000 J1820 
00001930 
00001640 
00001850 
OOOOlsbO 
00001870 
000 0 I860 
00001390 
OOoO 1900 
00001910 
00001 9£0 
00001930 
0000 19U0 
00001960 
00001960 
00001 970 
0000  19IJO 
0000 1 990 
00002000 
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CjM“mQN/9LGC<2/DISTA'N(75»75J»  I POLFL ( 7b , 75 1 
INTEGER  TEST 
C 

c 

c 

«Rl TE ( IQUT , 1 b5 ) NODES 
lbi  FORMAT (/*  SPATh,  NOOES=',lb) 

C 

C The  DISTANCE  MATRIX  OISrANU.J)  MAS  SEEM  CRtATED  I N SU5R.  SETUP. 

c 

C FLAG'S  ALGORITHM.  K is  THE  ITERATION  NUF'bER. 

c 

00  50  **1 « NODES 
00  50  1=1, NODES 
IF(I.ta.A)  GO  TO  50 

OISTIK  s OISTANCI/K) 

IF(OISriR.EU.9999)GO  to  50 
00  49  J=l, NODES 
IF ( I . EQ . J 1 G J TO  49 
IF(k.£U.J)GO  TO  49 
DISTkJ  s DISTAN(K.J) 

IFIDISTKJ.EU.INF1 Y)GO  TO  49 
TEST  = OISUKyOISTAJ 

c 

c 

IF (DI S T Am (I , J ) . L£ • T£S  r ) GO  TO  49 
OISTAN(UJ)  = TEST 
IPOLFL ( I « J 1 = IPQLFL ( I , K ] 

49  CON T I NuE 

50  CONTINUE 
RETURN 
End 

subroutine  rank(ngo,n,r,8,iii 

C R - ARRAY  TO  oE  RANKED. 

C ONLT  TMAT  PART  OF  R From  R(NGO)  TO  R(NJ  IS  RACKED 
C NGJ  = STARTING  3'UdSCR  I P T (IF  ARRAr  R 
C N = E.OING  SuHSCRIRT  OF  array  R 
c 9 = RAi.rEU  array. 

c ARRAY  II  KEEPS  track  OF  TmE  SUBSCRIPTS  OF  R 

C 11(11  SMOuLO  BE  INITIALIZED  Im  calling  PR0GRA4  isITM  11(11=1,  1=1. 2, 
C PRIOR  TO  FIRST  CAuL  OF  RANK. 


00002010 

00002o20 

00002030 

00002040 

0000205 0 
00002000 
000020/0 
00002090 
00002090 
00002 100 
0000211 0 
00002120 
00002130 
30002140 
00002150 
00002 1 bO 
00002170 
00002190 
00002190 
00002200 
00002210 
00002220 
00002230 
00O022U0 
00002250 
000022b0 
00002270 
00002290 
00002290 
00002300 
00002310 
00002320 
00002330 
00002340 
00002350 
000023b0 
00002370 
00002330 
00002390 
..03002400 
00002410 


C RE"3yE  C IN  CCJLJM-m  1 DF  STATE  RENTS  RERTaInI.mG  TO  II  1m  URCER  To  IN0UR00002423 
C II.  00002430 


Integer  r,s 

DIMENSION  R(1),B(U 
C OI.kENSIOim  II  Cl) 

DIMENSION  IKU 
NGl  s NGU+l 
NM  s NGU-l 
BINGO)  = R(NGCl) 


00002440 

00002450 

0'J0o24b0 

00002470 

00002490 

00002U90 

00002500 
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II(NGO)  = NG3 
*snGU 

OG  70  l*NGl,N 
IF(BCK)-RCI))  10,10,30 
10  3(1)  s w { I ) 

C 1 1 C i J = I 

II(I)  = I 

30  h = i 

GU  TO  70 

30  0U  50  U SNGO , A 
J = I-L  »NM 

IF  (BU)-R(  I ) ) bQ»e0,u0 
uo  d( J*i ) s BCJ) 

C IICJM)  S II(J) 

IICJM)  = II(J) 

50  CONTINUE 
BCJ)  s R(I) 

C II(J)  = I 

1UJ)  = I 
GO  TO  20 

to  a c j ♦ i ; = wen 

C llUtl)  = i 

lUJtl)  = I 

KSI 

70  continue 

RETURN 

END 

SUBROUTINE  CNTLRDCIUST) 


this  SUBROUTINE  READS  THE  CONTROL  FIL£  To  DETERMINE  ahICH  SHORTEST 
PATHS  SHOULD  BE  FOUND. 

VARIABLES: 

NETNUM(I)  : NUMBER  OF  I-TH  nETnQHK 

I Flag  ( I ) s 0 : FIND  all  ShURTEST  paths  for  „et,<QRk  I 

= 1 : F I NO  ALL  paths  FRGh  All  OTHER  NODES  TO  IDEST(I) 

= 2 J FIND  All  PATrtS  FROM  IOHIG(I)  to  all  other  NODES 

= 3 : FIND  SHURtEST  PATH  FROv  lQRIG(I)  TO  IOEST(I) 

IORIG(I)  : ORIGIN  OR  SOURCE  vOOE,  NETaURk  I 
IDEST(l)  S DESTINATION  NODE  FOR  vETaORK  I 
NETS(I)  : array  to  store  values  of  neinum 
ICNTLS  : VALUES  uf  IFLAG 

ISOURC  : VALUES  OF  IORIG 

I S INa  : VALUES  uf  idest 


lqhmon/io/in, iuut 

COMMON/BLOCKl/NlNtlb, INF  IN, K IMF  IN,  I T «U 1 6 , N£ TS ( «00 ) 
DIMENSION  Nt T NUM ( 7 ) 

ILAST  * 0 

30  READCIN,«0,END=9999) (NE  TftU  *(  I ) , IFLAG,  IuhIG,  IDEST,  1*1  ,7) 
aRITECI UU T, 40)( NET  NUN (l), I FLAG, IORIG, IUE3T,I  = 1,7> 


00003310 
00003530 
00003530 
00003540 
00003550 
0O0035eu 
000035 7 0 
0 0 0 0 3 3 3 0 
000035RQ 
00003600 
00003610 
00003630 
00003630 
00003b40 
00003630 
00003660 
00003e7o 
00003680 
00003690 
00003700 
00003710 
00003730 
00003730 
00003740 
00003750 
00003760 
00003770 
00003780 
0O003790 
00003800 
00003310 
00003330 
00003830 
00003340 
0 0 0 0 3 3 5 0 
00003860 
00003370 
00003830 
00003890 
O00039U0 
00003910 
00003930 
00003950 
00003940 
00003950 
0O0O3960 
0 000  3970 
00003990 
000O39QO 
00003000 
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40  F3RmaT(7(I3»I1»2I3)) 

00  50  1 = 1 , 7 

IF(NETNUMCI).lE.O)  GO  70  9999 

ilast*ilast*i 

50  NETSCILAST)  = NETNuM(I) 

GO  TO  30 
9999  REMIND  In 
RETURN 
6 NO 

SUBROUTINE  PR'IOUTIDISTAN, NODES,  I A x J 
CQMMQn/IO/IN.IOUT 

C0i4M0N/BLUCKl/NIriEt6,INFIN,MNFlN,  ITnOIo,  NETSC400) 

COMMON/ oL0C.<2/ 10  TAN (75, 75)  ,F  I LET  (75. 75)  , NO  TWA.,  (75)  , I NOE  XT  (75) 
INTEGER  UlSTAN ( 75, 75) » hIARCH(5) 

PRINT  OUT  FUR  FLOYO'S  ALGORITHM 

00  1 1=1,5 

1 HIARCH(I)=UNMASK(NETS(IAM),I) 

NETIO  = h I ARCH (5 ) 

•VfiITE(IOUr,lO)HlARCH 

10  format ( / ' MATRIX  FOR  xETrtORK  : ZONE  = ' , 1 5/ iux 'EOC  = 

1 , 30X 'GRUUP  =', I5/31X-SECTOR  = ' , 1 5/3 1 X 'UN  I I AREA  =',I5/) 

*G0=1 

23  *rtIT=*GU*14 

IFCKnlT.GT.  (ODES) an  I TsNOOES 

flRI T£ ( IOUT , 2b ) 

HR  I TECIOUT,  2b)  (NO  TRAN  (K  ) ,K=<GO,K’nIT) 

00  24  J=l, nodes 

..RI  Tt(IOUT,23)NOTRA;UJ),  (OISTANU.X)  ,x=KGU,R..IT  ) ,NOTRA(.(J) 

24  CONTINUE 

■>R  I TEtlUUT ,2b) (NOTRAN(K) ,r=RGO,RrIT ) 

X GO  = Kn I T ♦ 1 

IF  (KGO.LE.NOOES)GO  TO  23 
2b  FORMAT (bX, 1 5 1 tt ) 

28  F ORma T (2x14, 1518,14) 

RETURN 

ENU 

Subroutine  match(a,m,b,m,h, jj, jmatch) 

SUBR,  MATCH  COMPARES  ELEMENTS  OF  INTEGER  ARRAYS  A AND  8 FOR 
EQUALITY.  EQUAL  ELEMENTS  ARE  A MATCH.  RETURNED  VALUES  ARE: 
IMATCH  (=  0 = NO  match;  = l = A MATCH) 

II  4 JJ,  THE  I NO  ICE  S OF  THE  MATCHED  ELEMENTS  IN 


00003010 
00003020 
00003030 
00003040 
00003050 
0000  30b0 
00003070 
00003030 
00003090 
00003100 
000031 10 
00003120 
00003130 
00003140 
00003150 
00003 1 bO 
00003170 
00003130 
00003190 
00003200 
00003210 
00003220 
',1500003230 
00003240 
00003250 
000O  S2b0 
00003270 
00003230 
00003290 
00003300 
00005310 
00003320 
00003330 
00003340 
00003350 
00005360 
00003570 
00003380 
00003390 
00005400 
00003410 
00003420 
00003430 
00003440 
0000  5450 


INTEGER  a(m),  b(n) 
I MATCH® 1 
00  1 I = 1 , M 


A AND  8 RE  SPEC  T I VE  000  0 34  bO 
0 0003470 
0000  5480 
0000  5490 
00003500 
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_s_ 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


11*1  00003510 

OO  1 J=t,N  00003520 

JJsJ  00003530 

IF(A(1).EQ.8U))RETURN  00003540 

l CONTINUE  00003550 

JMATCM  s 0 00003560 

RETURN  00005570 

END  00005530 

FUNCrtUN  STaR£(vALUE.«GRQ,  I ) 00003590 

IMPLICIT  INTEGER  C A-Zo  00003600 

PlHtNSIJN  4ASa(5)»  P0.'fcR(5)  00003ol0 

LOGICAL  S ( R 1 P ( 5 ) , P , L *QRO  00003620 

DATA  FIRST/O/  00003650 

IF(FIRST.EU.1)G0  TO  10  00003640 

00003650 

CREATE  4 A 3 A S (ALL  3IT3  ON)  IN  5 OIFFERENT  Q I T FIELDS  TO  3E  LSED  FQR  00003660 

storing  or  extracting  corresponding  dir  fulls  in  a ..jrd.  beginning  oooo367o 

rtITH  TnE  HIGHEST  ORDER  3ITS  Tm£  NU«5ER  OF  dITS  IN  EACH  FIELD  ARE:  uOOOSbRO 

5 , 6,  6,  6,  9.  THESE  FIElOS  CORRESPOND  TO  VALUES  IN  THE  AkRAT  m I AhCh 0 0 0 0 36 90 


<hh ICh  are  zone  I.D.,  EOC  I.D.,  GROUP  I.D.,  SECTOP  I.D.  AND  nET*ORa 
(OR,  UNIT  AREA)  1,0.  These  5 VALUES  are  STOhED/ExIHACTED  IN  aOWD. 

2 first  = 1 

PQaER(I)  s 2*. 27 
PUaER(2)  S 2**21 
PO«ER ( 3 ) = 2**15 
PGaER ( 4 ) s 2**6 
POaER (5 ) * 1 
TnOS  * 2»*5-l 
4ASk(1)  * T n05 *PO«ER ( 1 ) 

TaUb  * 2**6  - 1 
MASK(2)  = T nOb*POa£R ( 2 ) 
m A 3 a ( 5 1 * TrtU6*PuwER( 51 
«ASa(U)  = T aU6  * POn£R  C 4 ) 
mASa (5)  s 2*  *9-1 
OO  1 J= 1 , 5 
1 STRIP(J)  = -aasrcJ) 

STORE  VALUE  IV  The  I-Th  dIT  FIELD  OF  aQRD. 


00003700 

00003710 

00003720 

00003730 

00003740 

00003750 

00<J037o0 

00003770 

00003760 

00003790 

00003800 

00003810 

00003820 

00003630 

00003840 

00003850 

00003860 

00003870 

00003880 

00003890 

00003900 


10  LhURD  * aORD 
Ps.NUr.STRIP(I) 

LAORO  = LrtORD.ANU.P 
A*VALUE*PU«ER(I) 

P s K 

P a LaORD.ANO.P 
STORE  * P 
RETURN 

Entry  unmaskuuro,  t ) 

LflORO».»UHD 


00003910 
00003920 
00003430 
oOOO  39u0 
00003950 
00003460 
0000  5970 
00003980 
00O03990 
00004000 
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P S IrtUWO.A'(D.STRIPCI) 
UNMASK  = H 

UNMASK  = UNKASK/POk£R(I) 

RETURN 

END 


U 0 0 0 A (J 1 0 

oooouoau 

00004050 
0 0 0 0 « 0 a 0 
00004050 


subroutine  4li.nET 

REMOVE  Th£  ABOVE  *C'  TO  USE  &LLNE1 
also  REPLACE  'STJP'  AlTH  'RETURN'. 


00000010 

TO  USE  ALLxET  as  A SUBPOjTIMt  CA.LED  BY  N£  T aORk  . 0 00  0 0 0 2 0 


ALSO  REPLACE  'STOP'  AlTH  'RETURN'.  00000030 

00000040 

00000050 

RT I. Caj.PO«95b.jH0.ALLNET. FORT  OOOOOOsO 

00000070 

THE  OBJECT  OF  THIS  STEP  IS  TO  CREATE  T h£  INITIAL  DISTANCE  MATRIX  OOOOOOBO 
OF  DIPECTLT  LINKED  UNIT  AREAS  AS  2 MATRICES  STORED  BY  ROaS  ON  OUTPUT  OOOOOORO 
DEVICES  91  AND  92.  Th£  MATRICES  AWE  INPUT  TO  Th £ NEXT  STEP , DA*A Ge  3000U1J0 
ASSESSMENT,  «*HlCH  ADJUSTS  ThE  DISTANCES  PEW  'BOS*  CODE  VALUES  And  00000110 
ADDS  THE  CORRESPONDING  ELEMENTS  IN  EACH  ROa.  00000120 

file  53  CONTAINS  RESULTS  OF  SHORTEST  Paths  COMPUTED  UITHIN  Each  00000130 
UNIT  AREA,  I.E.  NETWORK,  ANO  The  IDENTIFICATION  OF  ThE  NODES  IN  EACH  00000140 
vET.tORK.  00000150 

TaO  UNIT  AREAS  A AND  8 ARE  LINKED  IF  THEY  RAVE  ONE  OR  «0r£  BOuNOAH 000 0 0 1 60 
NODES  In  COMMON.  00000170 

LINKAGE  OF  A ANO  d IS  DETERMINED  BY  COMPARISON  OF  THE  LIST  OF  00000130 

NODES  OF  EACH  NETaORK.  EQUAL  NODE  I.D.'S  UCCJH  FOR  COMMON  NUDES  00000190 

UN  Th£  BUNCAWy  BETaEEN  THE  T.vQ  NET«OWK3.  00000200 

IF  NETaORKS  A AND  B ARE  CIRECTOT  OINKED,  ThE  DISTANCE  FWQm  A TO  B 00000210 
IS  DEFINED  AS  THE  AVERAGE  OVER  SHORTEST  PATHS  FROM  ALL  NUDES  IN  A TO  00000220 
THt  COMMON  dOUNDARY  NQOtS  PLUS  Th£  AVERAGE  OVER  SHORTEST  PATrS  FRO"  00000230 
The  COMMON  BOUNDARY  NODES  TO  ALL  NODES  IN  6 , 000002U0 

00000250 

DIMENSION  TS ( 400 ), NOO I NA( 75), MOO  I * 8 (75). SPA  THAI  75, 75), SPAT *6(75,00000260 

1 75)  00000270 

INTEGER  A2 3llUN ( 40 0 ) ,B2BOUN(«00) ,BOUN2A(uuO) ,BOLN2b(40  0) , 000002BO 

1 T 080UN/9 1 / , FRMBUO/92/,  SP A Tha , SP A ThB  00000290 

COMMON/ SHIFT/ 1 r.,016, INF  IN, A2BUUN,028OUn,8OJN2A,6OUn29  00000  30  0 

00000310 

IT.vOlo  = 2**16  00000320 

K INF  I N s 9999  00000330 

INFI„  s KINFIN»IT«Ol6  00000340 

REadC3u)maxnET,NETS  00000350 

RE  A I NO  30  OOOOOSoO 

NUMNET  s MAX.NET  00000  3 70 

00000  3B0 

I/O  DEVICES:  00000390 

180  s 60  00000400 


DO  10  Ns  1,40  0 
A2B0JN(N)sINF IN 
60UN26 ( N ) s HF IN 
B2B0JN(N)SINFIN 
10  dOUN2A ( n ) s INFIN 
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INITIALIZE  FILES  70  At.O  80;  ASSUMPTION  IS  THAT 

distances  initially  egual  Infinity  (=oppp). 


ALL  DlKECTLY  LI  (*£00000051 0 
00000520 


00000530 
OOOOObAO 
00000550 
00000560 
OOCOOS70 
00000580 
00000500 
OOOOOoJO 
0 00  0 Oo 1 0 
00000620 
OOOOOoiO 
00000600 
00000650 
00000660 
00000670 
0000o660 
OOOOObPO 
00000700 

■ 00000710 

P£AO  FILE  53  FQR  the  results  3F  ThE  SnUKTEST  PAThS  COMPUTATION  FJP00J0O72) 
UNIT  APE  A A.  THE  FIkST  8ECORO  CONTAINS  -.0.  OF  NETaQPKS,  AND  NET  I D * S 0 00 0 J 7 3 d 

000  00  7 a 0 

PE  AO ( 5 3 ) InOExa,  »R\iaOA,iNOOlNA,SPATHA 
-ETnUA  — UN-ASKd  <C£«A,5) 

A28UUN ( 1 A J s NET  40A 
B0UN2S  ( I * ) * <£T.OA 
M A X A J «CDJnA(  .a-iODAJ 

P«  I'.T  1 xH2.  I a,  .£  T .l  A.  NR.UOA,  ( 'lOOINA  CL  N)  , LH=1  ,NWi*UQA  ) 

1*P2  Fuh^aT  ( 30iA.il  , 5 « j . I 4 }) 


DU  11  4S1,U00 

*PITE(I70)N* CINFIn, 1 = 1 ,U00) 

11  «HITECI80)N, (INF  In, 1 = 1, u oo) 

PEaIND  170 
HE.-iIND  180 

initialize  loop  i jdex  to  count  unit  area  *a' 

IA  s 0 
1 IA  s IA*i 

ifua.gt.haxnetj  go  ru  200 

INITIALIZE  the  yECTOH  OF  DIPECTlY  LINKED  DISTANCES  fop  NET.-.OPK  A 
the  vECTUP  IS  The  IA-TH.  RUa  of  THE  DISTANCE  ha  T P I X 

RE AO (1 70, ENOS 1 00) INDEX A, A290UN 
PEADII80)  I NDE  X A , dOuNZB 


Nt  T rtJM*  a » \ v.fc.  i .^*HAA£D  «ITm  nETaOPkS  8 (5rIAtl,lAA2. 
...*  <t  T«  5 . 

1 8<  I A 

2 lu*Io»l 

IF (IB.GT . HAXNt TJGU  TO  100 
PtA0(I70»tND*l00Tl  .L'f  «B  , d2BUU*l 
0EA0(I8oHnJEx8,b3un2a 


00000751 
00000780 
00000770 
000007*0 
000OO7PO 
00000500 
00000810 
00000520 
0 J 0 0 0 A 3 u 
OuuOOhuO 
000008 jT 
00000860 
0000057 0 
00000850 
500G05PO 
OOOOONOO 


OOOOOOIO 

A3Sd“E  initially  T~aT  A AND  d APE  NOT  LINKED  AND  Th£  DISTANCE  A TU  rf  0)0 0)P2> 


IS  INFINITE. 

B2HUUN 1 1 A ) s INF  IN 
80UN2ACIA)  s IuF IN 

P£Au  Th£  ShUPTEST  PATHS  PESULTS  FUH  n£T«UHK  8 
PE AO ( 5 3 > ENO* 100)  I NOE  XB # NRNUOB » NO0 1 Nb » SP  A f h8 
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OOOOOP 30 
OOOOoPuO 
OOOOOP8C 
OOOOOPbO 
0 0 0 0 0 P 7 0 
00000080 
OOOOOP  )0 
00001000 


I 

I 

I 


I 


I 

I 


NETNOB  2 Ui*4ASN(IN0EX8,5) 

82B0JNUB)  = NETNUS 
80UN2A(I8)  r NETNOB 

C PH  I. 'll  149a,  I8,NET«0B,NRN00B,  (NUDINB  ( LM)  , L ”2  1 , NR.NODB  ) 

c 

c l *ill  be  used  to  cuunt  the  number  of  ha rcHto  nodes  in  a and  a. 
c 

L20 

A3BOU  2 Q. 

B2SUU  2 0. 

80U28  2 Q. 

80U2A  = 0, 

JGU=1 

MAX3  2 NUUlN8(NHNOD8) 

DO  4 1 3 1 , nW  NOD A 
ITHNOD  3 .MOO  I.M  A C I J 
C PRINT  1492, l.ITHNOO 

IF (ITHNOO.lt. NODINB(JGO))GJ  TO  4 
DO  J J2JG0, NRNOOB 
JTHNOO  3 NOUINB(J) 

. IF  (ITHNOD. GT.MAXB. OH. JTHNOO. GT . VIA < 4 ) GOTO  45 

C PRINT  149J, J, JTHNOO 

1495  FOR*Ar (5X31 5,  ' J,  JTHNOO') 

IF (JThNuO.GT.ITHNOC ) JGO=J 
IFUThNOO.GT.IThnJDJGO  TO  4 
IF ( J ThnOO .LT . I ThNOO ) GO  TO  5 
C 

C A HATCH  BETWEEN  ITH  NODE  IN  A AND  J-TH  nQDE  1 9 B 

c 

L 2 u»l 

c the  following  give  o i 3 t a to  b as  ( a ro  bound.  » bound.  tu  bj. 

A3B0J  2 A2BUU*C0LSUM(SPATHA, NRNOOA, I ) 
oOuBb  3 80U33  ♦ ROwSUMCSPATHfl.NHNOOB, J) 

C The  FOLLOWING  GIVE  DIST  8 TO  a as  (8  TO  SOUND.  ♦ BOUND . TU  A). 
o360U  = B3B0U  ♦ COLSUMCSPAThB, nRhODB*  <!  1 
B0U3A  2 SOU  2 A . Ro  'iSUH  ( SPA  T H A , nRI.ODA  , 1 ) 

JGU  * JGU»1 
GO  TO  4 
5 CUNTINUE 
4 CONTINUE 

C L 2 0 s NO  NODES  IN  CUMmqn  bETwEEN  NETnuRkS  A AND  8. 

45  IF(L.£Q.O)GO  TO  4b 

IA380U  2 A3BOJ/(L»nhnOOA)  ♦ 0.5 
A3B0UN  (18)  2 IA3bOU»I  T.'.Ulb  ♦ NET  NOS 
IB0U38  2 BUU3B/(L*NRNODs)  ♦ 0.5 
BOUN36  (18)  2 I fl(jU33  » I T*0  lo  ♦ NETNOB 
C 

C THIS  COMPLETES  Th£  IB-Th,  HO*  OF  ARRAY  'TO  OUu'i.'  iP  1 
c Element: 
c 


OoCOlOl  ) 
00001030 
00001 050 
00001U40 
00001050 
00001 060 

0 0 0 U 1 o 7 0 
00001080 
00001 090 
0000 1 100 

00001  1 1 0 
0000  1 1 30 
00001 1 50 
00001 1 4U 

00001150 
00001 1 bO 
0000  1 1 70 
00001 l 30 
00001 190 
00001 200 
0000  121  0 
00001220 
0000  1230 
00001240 
00001250 
0000 12S0 
0000 1270 
00001230 
00001290 
0000  1 300 
0000 1 510 
00001320 
00001330 
000013-0 
00)0 1 550 
0 J 0 0 1 is 
00  0 ) 1 3 • . 


I 

I 
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AO-A036  265 


UNCLASSIFIED 


RESEARCH  TRIANGLE  INST  REStARCH  TRIANGLE  PARK  N C OPE— ETC  F/G  15/3 
LOCAL  EMERGENCY  OPERATING  SYSTEM  - LEMOS. (U) 

JUL  76  J H DUNN.  R N HENOKY.  R 0 LYQAY  DAHC20-73-C-0253 


RTI-AAO-873 


NL 


END 


DATE 

FILMED 


*5^7? 


I 

I 

I 

I 


I B2B0U  x B2BJU/IL<«NRN0DBJ  t 0.5 
B280UN(  14)  x !S2dUU*irn01b  ♦ NETnOA 
C 

C THIS  COMPLETES  the  18-Tn.  SO*  OF  array  'FROM  dOUNO.'  UP  TO  The  I a Th . 
C Element: 

c 

IBOU2A  s 80u2a/(L»NRn00A)  t 0.5 
30UN2A ( I A } s 1 BOU2  A • I T wO lb  ♦ \ETNU* 

4b  IN0EX8  3 NE T NU8  » l TwO  l b ♦ Id 

wRITEd 71 ) IN0Exa,B28UUN 
oRiTEdan  tNOExa.aouNEA 

c 

C CONTINUE  processing  unit  AREA  a against  The  next  unit  arFA  d. 

c 

GO  TO  2 

C 

C 100  « HAvt  PROCESSED  all  8 NETWORKS  ‘GAINST  NfiTwOfiK  A. 

C WRITE  ThE  RESULTS,  A2BUUN ( J ) A <0  60UN28 (J J , J3 1 , 2 . . , , u n£ T wORk S , 

C UN  FILES  'T080UN'  ANO  'FRM80U*. 

c 

100  REWIND  53 
REwINO  70 
REwINO  80 
REWIND  71 
REwINO  SI 

INOEXA  s N£  TNJA* I T *0 1 6 t IA 
WRITE! T080UN)  INDEX A , A2dOU  N 
WRITE (FRMBOU ) INDEX A, 8UUN28 
C 

C advance  file  53  TO  THE  NEXT  UNIT  area  A AHlCh  is  TO  BE  PROCESSED 
C AGAINST  all  other  a UNIT  AREAS  EXCEPT  8<A, 

C 

DU  12  L*1 » I A 
12  REAO(53,ENO*200J 
C 

C PART  OF  The  NEXT  VECTOR  A has  ALREAOr  BEEN  PROCESSED  AS  A vECTUR 
C 8 UP  TO  The  IA-Th.  ELEMENT  AND  The  RESULTS  /.£R£  wRITTEn  UN  FILES  171 
C ANO  1 8 1 

C SO  READ  THOSE  FILES  FOR  VECTOR  A INSTEAD  OF  READING  1 70  ANO  180. 

C TO  DO  THIS,  FlIP-FLUP  ThE  VALUES  FUR  180  AND  170,  181  AND  171. 

170  3 141-170 

171  3 i«i-I71 

180  3 161-180 

181  3 lb  1 -1 8 1 

GO  TO  1 

C 

C 200  3 ALL  unit  area  NETWORKS  HAVE  been  PROCESSED  AS  tET-ORk  'a* 

C FOR  'A'  3 i,  2,  ...  , I A— l . 

c 

200  REwINO  53 


00001510 
00001520 
00001530 
00001540 
00001550 
0000l5b0 
000O 1 57o 
00001580 
00001500 
0000 1 bO  0 
OOOOlblO 
00001620 
OOOOlblO 
00001640 
90001650 
0000 1 bbO 
00001670 
0000 1 680 
00001 690 
00001700 
00001710 
00001720 
00001730 
00001740 
00001750 
00001760 
00001770 
00001780 
00001790 
00001800 
00001810 
00001820 
00001830 
00001840 
00001850 
00001860 
00001870 
0 0 0 0 I H 8 0 
tlOOOlSRO 
■10001900 
00001910 
00001920 
00401930 
00001940 
00001950 
00001960 
00001970 
00001980 
00001900 
00002000 
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c 


REnlND 

70 

00002010 

RE" INO 

80 

00002020 

RE«INO 

71 

00002030 

REoINO 

81 

00002090 

R£"IN0 

TOBOUN 

00002050 

REnINO 

FRrtBOu 

00002060 

00002070 

The  INITIAL  0I3TANCE  matrix  OF  OIRECTLY  linked  nodes,  E.G.  unit  00002080 
AREAS,  NO"  EXISTS  AS  2 FILES  OF  R0"S  ON  OEVlCES  91  AND  92.  Tm£  F ILES00002090 
ARE  DISTANCES  from  all  NODES  IN  A TO  ThE  boundary  bETwEEN  A a NO  8 00002100 
ANO  TmE  DISTANCES  FRO*  THE  BOUNDARY  BETnEEN  A ANO  B TO  ALL  NODES  IN  B00002110 


USING  ARRAY  B0UN2A  FUR  "ORKING  STORAGE  IN  SUBR.  NETPRN; 

CALL  NETPRN{MAXNET,BOUn2A, TOBOUN, l) 
call  NETPRN(.MAXNET,B0UN2A,FRMB0U, 1) 
call  N£TPRN(MAXNET,BOUN2A,TOBOUN,2J 
call  NETPRN(MA*n£T,80UN2A,FRMBOU,2) 

STOP 

TO  USE  ALLNET  as  A SUBROUTINE,  REMOVE  'STUP*,  AND  REMOVE  'C'  IN 
•RETURN*  CAROl 
RETURN 

EN0 

function  ru«su'i(distan,nooes,ro*num) 

THESE  2 FUNCTIONS  ARE  APPLIED  TO  TH£  DISTANCE  “ATRICES  CUMPUTED  FOR 
THE  INDIVIDUAL  UNIT  AREA  n£T*ORKS. 


COMPUTE  ThE  SU9  OF  DISTANCES  FROM  NUDE  ROwNUM  TO  A».L  NOOES. 

INTEGER  OISTAN(75,75),RO"NU.a,COLNUM 
RO"SU*  * 0. 

DO  1 J*l, NODES 

1 RO"SUM  3 RIJ«3UM  ♦ 0ISTAN(H0"NU«, J) 

RETURN 

ENTRY  COLSUHIDISTAn,  NU0E3,  COL'YUN) 

COMPUTE  THE  SUM  OF  DISTANCES  FROM  ALL  NODES  TO  NOOE  COLNu<*. 

COLSUM  s 0. 

DO  2 1*1. NODES 

2 COLSUM  3 COLSUM  ♦ OI3TAN(I,COLNUM) 

RETURN 

EnO 

SUBROUTINE  1ETPRN(N,LINEIN,I0,I"HICH) 


00002120 
00002130 
00002 1 RO 
00002150 
00002160 
00002170 
00002180 
00002190 
00002200 
00002210 
00002220 
00002230 
00002290 
00002250 
00002260 
00002270 
00002280 
00002290 
00002300 
00002310 
00002320 
00002330 
OOOU23UO 
00002350 
00002360 
00002370 
00002380 
00002390 
00002900 
00002R10 
00002920 
00002930 
00002990 
00002R50 
30002960 
00002970 
00002980 
00002990 
00002SO0 
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I 


I 

I 

I 

f 


It 


C I/IMICM  « 1 « P R I ;<•  T OUT  DISTANCE  VALUES  - UPPER  mALF  C F rtOMD 
C IaHICm  S 2 3 PRINT  OUT  ROLlCV  VALUES  - LO»EK  HALF  OF  *0»D 

c 

DIMENSION  lI'.EI  n(  1)  .LINOUT  (31) 

CUHMON/ShIFT/ITa016, IwFIh 

c 

C OUTPUT  DEVICE  TOUT: 

I OUT  a } 

JGO*l 

5 »»lTE(luur,3) 

3 FCHMAT(lwl) 

t «EAO(IO,ENO*100)INOEx,  (LINE  IN  CU  .L»l.N) 

JArtlT  a JGOAiO 

IF(JK<*lT.GT.N)jK*ITaN 

Lt*0 

00  a L*JGO , N 
EL  a LL*1 

LlNOUT(LL)»lI NEIN(L) /I  T.iOlo 

IF(I.VMlCM,£a.2)Ll^0UT(LL)3LI'VEI'»(L)  - LINUU  T ( LI)  * i T*U  1 6 
tt  CONTINUE 

1 NOE XI  a UN, MASK!  INDEX.  5) 
aPITCdUUT.2)  IN0EX1,  (LlNUUT(L).tai  ,LU 

2 FUHPATUXI3.31U) 

go  ru  i 

100  REMIND  10 

JGO  a JKaITM 
IFUGO.GT.NJRETUHN 
GO  TO  5 

End 


00002510 
00002520 
00002530 
000025u0 
00002550 
00002560 
00002570 
00002580 
00002500 
000026u0 
00002610 
00002620 
00002630 
000026U0 
00002650 
00002660 
00002670 
00002680 
00002600 
00002700 
00002710 
00002720 
00002730 
00002  7u0 
00002750 
00002760 
00002770 
00002780 
00002790 


r 


r 
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c 

C TO  USE  'DAMAGE'  AS  A SUBROUTINE#  REMOVE  'C'  IN  SUBROUTINE  CARO  ANO 
C CHANGE  'STOP'  TO  'RETURN'. 

c 

C SUBROUTINE  damage 
C RTI.CAI  .P0«95b.J-»U. DAMAGE. FORT 
C PERFORMS  DAMAGE  ASSESSMENT  OF  UNIT  AREAS. 

C ASSIGN  Ntn  DISTANCES  USING  BOS  VALUE  DN  RESOURCE  FILE 
C THE  RESOURCE  FILE  IS  ON  FT61F001. 

c 

IMPLICIT  INTEGER  U-Z) 

DIMENSION  NETS l 400 J, DAMAGE (400), A2BOUNC 400) , BOUNDS (400), 

1 D1STAN(400) 

COMmOn/SH IF  T/ 1 TnO lb, INFIN 

l T aOI b s 2 • » lb 

KlNFINa9999 

INFIN  « A INF  IN* I T aOI b 

REAO( SO JMAxnET » N£ TS 

HEaINO  SO 

NUMNET  * MAXNET 

00  to  I*t, NUMNET 

NETNO  a UNMASK(NETS(I)>S) 

IFCNETNO.LE.OtGO  TO  20 
8 REAO(bl,ll,ENO*20)  ZONE  ID,  UN  l T 10, LUC , BOS 
U FORMAT (2IJ, 12,SXI1) 

IFCLUC.NE.5tG0  TO  8 
C 

C FINO  A match  BETWEEN  the  UNIT  AREA  ID.  UNlTID,  AND  SOME  nETauRa, 

C NETNO.  IN  THE  LINKS  file. 

c 

IF ( Nb  TNO.LT ,UNI T ID ) GU  TO  8 

c 

C redundancy;  ASSUMED  sort  on  NETNO  ANO  JNITIO  should  PRtCLUOE 

C NET  xQ>UN I T 10 

c 

IF CNE TNO.GT .UNI T IDtGO  TO  Iv 
OAMAGECn  a 2*  * ( oOS- 1 ) 

10  CONTINUE 

C 

C RE  AO  A2S0UN  and  B0Uw29  FILES  AND  COMPUTE 
C 0 ( I*J)*A2bUUN(J )*OAMAGE(I)  * B0UNZ8(J)*DAMAGE(J) 

C 

20  REaInO  bl 

DO  55  I«l. NUMNET 

RE AO (9 | ) INOEXA, A280UN 

RE AO (92 ) INDEXA.BOunZB 

DO  54  J«l. NUMNET 

IA260U  a A28 JUN ( J ) / 1 T aO 1 a 

DISTANCJ1  a INFIN 

lFCIA2BOU.EO.KlNFlNjr,0  TO  54 


00000010 
00000020 
00000030 
00000040 
00000050 
OOOOOObO 
00000070 
00000080 
0 00  0 u 090 
000001  JO 
J00001  10 
00000120 
00000130 
0000014U 
00000150 
OOOOOlbO 
00000170 
00000180 
00000190 
00000200 
00000210 
00000220 
00000230 
00000240 
00000250 
OOOOOZaO 
00003270 
00000280 
00000290 
00000300 
00000310 
00000320 
00000330 
00000340 
00000350 
OOOOQ3bO 
00000370 
00000380 
00000390 
00000400 
00000410 
00000420 
00000430 
00000440 
00000450 
000004b0 
00000470 
00900480 
90000490 
00000500 
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I 

I 


I 

I 


I80U2B  s BUUN2B (JJ/IT*016 
POLIJ  * UNMA3K (NETS ( J)  ,5) 

DISTAN(J)  a UA2BOU*OAMAGE(I)  ♦ 1 90U2B ‘DAMAGE IJ))*ITaU16*POLIJ 

54  CONTINUt 

55  *fiITE(9i) INOtXA,Ul3TAM 
REMIND  91 

REmIM)  92 
MEmInO  95 

CALL  NfcTPHN ( 4AANET  > A^dOUN* 95* 1 ) 

CALL  NETPHN(MA*\iET,A290UN,9S»2) 

STOP 

RETURN 

END 


00000510 

00000520 

00000550 

00000540 

00000550 

00000560 

00000570 

00000560 

00000590 

00000600 

00000610 

00000620 

00000650 

00000640 

00000650 


? 


i V 

% I 

i 


r 


\ y 

I 


K 


l 

r 
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C RTI.CO3.PUO956.JaD.Fl0Y0.F0HT 
C 

C TU  USE  'FLOYD*  AS  A SUBROUTINE,  BEHOVE  C IN  SUBROUTINE  CARO  AnO 
C CHANCE  'STOP*  TO  •RETURN* 

C 

C SUBROUTINE  FLOYO 

C 

C PROGKAM  FUUYO.  COMPUTES  SHORTEST  paths  in  OIST(I,J)  USING 
C FLOTO'S  ALGORITHM  AOAPTEO  TO  aQHK  aITh  ONLY  2 RO*S  Of  0 
C IN  COHE  AT  UNE  TIME.  0(1, J)  IS  STORED  dY  ru.*S  ON  CEvICE  93. 

C THE  SOLUTION  MATRIX  IS  nRITTEN  TO  DEVICE  55. 

C 

INTEGER  U3InG(O00),CHANGE(O00),TE3T,013TIk,DISTIJ,0ISTKJ 

IT *01 6 • 2**16 

K INF  IN  a 9999 

INflN  a K I NF IN*IT*016 

ITO  « 70 

ISO  > SO 

NEXUSE  > 60 

KGOat 

REAO  (30)  NUMNET.NETS 
REnINO  30 

39  00  12  Ia|,NUMNET 

R£A0(93)1N0EXI, CHANGE 
12  aRITE(170)1V0EaI, Change 
REaINO  93 
RLNlNO  170 
C 

C INITIALIZE  the  FIRST  'using*  ROh 

c 

REA0(I70) INDEX  I , USING 
REhInU  170 

aRITE (NEXUSE) INDEX l, USING 
REaINO  NEXUSE 

c 

C BEGIN  FLOYD'S  ALGORITHM. 

c 

00  ttO  AaKGO, NUMNET 

c 

C EVALUATE  all  CURRENT  SHORTEST  PATHS  0( I, J)  USING  each  node  a 
C IN  TURN  AS  AN  INTERMEDIATE  NUDE  IN  ThE  PATH.  IF  THE  DISTANCE  USING 
C NODE  K IS  LESS  THAN  ThE  DISTANCE  mIThOUT,  INSERT  NUDE  * INTO  The 
C PATH.  0CI,J)  IS  STUHE0  BY  ROaS  Oi<  OEvICES  93  AND  170. 

c 

C 

C READ  In  THE  *U3INC«T0-CmAU6E ' »U»,  I.E.  The  *-Th  Rom 

c 

RE AO (NEXUSE) INDEX A, USING 
REaINO  NEXUSE 
DO  30  I a | , Nu  RNET 


00000010 

00000020 

00000030 

OOUOOOttO 

00000050 

00000060 

00000070 

00000060 

00000090 

00000100 

00000110 

00000120 

00000130 

OOOOOluO 

00000150 

00000160 

00000170 

00000180 

00000190 

00000200 

00000210 

00000220 

00000230 

OOOOOZttO 

00000250 

00000260 

00000270 

00000280 

00000290 

00000300 

00000310 

00000320 

00000330 

000003UO 

00000350 

00000360 

00000370 

00000380 

00000390 

OOOOOttOO 

OOOOOU10 

00000920 

00000U30 

OOOOOttttO 

00000050 

00000O60 

00000070 

00000O80 

00000090 

00000500 
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c 

C RE AO  IN  rut  'T3-8E-CHANGED*  RJ«,  I.E.  THE  I-Th.  RCU 

c 

REAO 1 1 70 ) INOEX I , CHANGE 
I7U.EU.DC0  TO  21 
C 

C OISTANCE  I TO  A IS  STORED  AS  AN  INTEGER  IN  UPPER  HALF  OF  THE  aOMO 
C Chance (k ) t 
c 

OISTIn  a Change (k)/ITh01o 
IF (DISTIk.E3.kI »F INIGO  TO  21 

c 

C UPOATE  The  J a 1,2. ....NUMNET  ELEMENTS  IN  THE  R[)w,  I,  bEiNG 
C CHANCED. 

C 

00  20  J*1 * NUMNET 
IFU.EQ.JIGO  TO  20 
IF(I.Ew.J)GO  TO  20 
OlSTKJ  a USINGCJl/lTAQlb 
IFIDISTkJ.E'I.KINFINIGO  TO  20 
TEST  a OISTIK  ♦ OISTKJ 

IF (TEST . LT. CHANGE! Jl/ITA0l6)CHANGE(J)aTEST*IT»0l6*(CHANGE(K) 

1 - OISTlKMTwOlb) 

C 

C STORE  THE  nE»  OISTANCE  IN  ThE  UPPER  HALF  OF  /. OPO  CHANGED)  Ano 
C STORE  THE  NE«  POUCT  IN  THE  lOhER  half, 

C 

20  CONTINUE 
C 

C "RITE  THE  CHANGED  ROn  ON  DEVICE  ISO; 

C 

21  wRITEUSOUNOCXI,  CHANGE 
C 

C SAVE  THE  SOLUTION  HATR1X  ON  FILE  55 
C 

IF (K. EG. NUMNET ) wRITE (55) DOEX I, Change 
IFII.nE.kIGO  TO  29 
00  26  lladNUMNET 
26  JSINGdll  a CHANGE  (It) 

29  IF ( I .NE .K*l )G0  TO  JO 
aRITE(NEXUSE) INDEX  I, CHANGE 
HEnInO  NEXUSE 

30  CONTINUE 

c 

C FLIP-FLOP  THt  I/O  DEVICES: 

c 

170  a 150  • 170 
160  a ISO  - 160 
REMIND  70 
HEwINO  60 


90000S10 
00000520 
OOOOOSJO 
90900540 
00000550 
OOOOUSbO 
OOOOOS  7 0 
00000560 
00000590 
00000600 
00000610 
00000620 
00000630 
00000640 
00000650 
UOOOOooO 
00000670 
00000660 
00000690 
00000700 
00000710 
00000720 
00000730 
00000740 
00000750 
00000760 
00000770 
00000790 
00000790 
00000600 
00000810 
00000820 
00000830 
00000840 
00000950 
00000860 
00000670 
00000680 
00000890 
00000990 
00000910 
00000920 
00000930 
00000940 
00000950 
00000960 
000009TO 
00000980 
00000990 
00001000 
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N 


I 


00001010 
00001020 
000010  JO 
00001040 
00001050 
00001060 
00001070 
00001060 
00001040 
00001100 
ooooi  no 
00001120 

I 

l 

f 


r 

1- 

r 


4 

r 


I 


rtT:.caj.pQ(*9b6.jwo.rvuw£c.Fo«r 

TO  use  TVLREC  AS  A SUBROUTINE,  REMOVE  'C*  IN  SUBROUTINE  Caro  and 
REPLACE  'STOP'  alTM  'RETURN  * , 

SUBROUTINE  TVLREC 

PROGRAM  TvLKECj  WRITES  RECORDS  FOR  TVLRtC  FIlE  ON  DEVICE  9«. 

implicit  integer  t a-zj 

DIMENSION  NETS (400), 01  ST AN (400) 

ITnOIo  * 2»»lb 

AS3UME0  VELOCITY  uF  TRAVEL,  SPEED,  IS  45  mIl£S/hUUR. 

SPEED  * 45 

READ! JOJNumnET, NETS 

REnINI)  30 

TRAVEL  TIMES  ARE  TRAVEL  TIMES  BETWEEN  UNIT  AREAS. 

THEREFORE  the  "LEvtL  of  DATA"  IS  5 CORRESPONDING  TO  LOWEST  LEVEL. 

LVLDAT  * 5 

00  50  1*1, NUMNE  T 

NETNUM  * UNMASK (NETS ( I ) ,5) 

30  RE AC(55HNDEXI,DIST AN 

ORIGIN  * unmask ( INGE* I , 5 ) 

IF(09IGIN.Nt.NETNuM)G0  TO  30 
TVLCQO  * 2 

REF  ZUN  a UNMA3K(NETS(I), 1) 

EJC  I.D.  13  UNMASK (NETS(I) ,2)  - NUT  REQUIRED  IN  TVLREC  RECORDS. 

RGHOUP  a unmask (NE TS( I ) , 3) 

RSECTR  a unMASk(nETS(I),4) 

kRI TE(94,9)REFZ0N, RCROUP,RSECTR,LvLDAT ,REFZ0n,0R IGIn,  T VLCUO 
9 F0RMAT(312, II, 4X, 213,11) 

TVLCUO  a 3 
00  49  J*l, NUMNE T 
IF ( I .EQ. J )G0  TO  49 
DI3T1J  a 0 1 S T An ( J ) / 1 T-<0 1 a 

0ISTANCE3  are  IN  MILE3M0,  I.E.  DISTIJ  * 2b  a 2.b  MILES,  FOR  EXAMPLE 
TVLTIM  a (DISTIJ*10)/SPEEQ 
ZONE  IS  ZUnE  OF  ThE  DESTINATION  AREA. 

ZONE  * UNMASK(NET3(J),1) 
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00000010 
00000020 
00000030 
00000040 
00000050 
OOOOOObO 
00000070 
00000080 
00000090 
00000100 
00000110 
00000120 
00000130 
00000140 
00000150 
OOOOOlbO 
00000173 
00000180 
00000190 
00000200 
00000210 
00 000220 
00000230 
00000240 
00000250 
00000250 
00000270 
00000290 
00000290 
00000500 
00000310 
00000320 
00000330 
00000  54(1 
OuOOu  530 
00000350 
00000570 
00000380 
OOOOu  540 
0C000U00 
00000410 
00000420 
,30000430 
00000440 
00000450 
0o00Q4b0 
00000470 
00000480 
00000490 

00000500 
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I 

I 

I 


C OEST  13  THE  IDENTIFICATION  OF  TnE  DESTINATION  AREA.  0U0G0510 

C 00000520 

DEST  * UNMASK (NETS (J ) ,5)  00000550 

aHITEIRR, 10 ) REF ZUn, RG»OUP, WSECTR, LvLDAT, TVLTIM, ZONE, OEST, TvLCOO  00000540 
UR  CONTINUE  00000550 

10  FUNMaT(J12,11,I«,2I3,H)  00000560 

50  CONTINUE  00000570 

STOP  00000580 

C 00000590 

C TO  USE  AS  A SUBROUTINE,  REPLACE  'STOP'  KlTN  'RETURN'  00000600 

C 00000610 

C RETURN  00000620 

C 00000630 

END  000006<i0 


I 


I 

I 

I 


